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INTRODUCTION 


This  final  report  documents  work  performed  by  Computational  niysics,  Inc.  (CPI)  for  the 
Space  Science  Division  of  the  Naval  Research  Laboratory  (NRL)  in  support  of  contract  numbtf 
N00014-88-C-2281.  Extensive  documentation  has  been  provided  throughout  the  contract  period 
in  the  form  of  nxmthly  reports,  as  well  as  more  formal  reports  detailing  the  completion  of  key 
research  tasks.  The  series  of  a{q)endices  contained  in  this  report  summarize  the  technical  work 
performed  in  key  research  areas.  In  addition,  the  following  section  gives  a  brief  summary  of  the 
general  scope  of  the  entire  contract 


OVERVIEW  AND  SUMMARY 

The  work  performed  on  this  contract  falls  into  the  following  general  categories: 

1 .  satellite  optical  instrumentation 

2.  analysis  of  optical  data  from  satellite  experiments 

3.  analysis  of  solar  spectra  and  optical  images  from  satellite  experiments  to  investigate  physical 
processes  taking  place  in  the  solar  atmosphere,  particularly  the  corona  and  transition  zone 
region 

4.  modeling  of  the  chemistry  and  dynamics  of  the  middle  atmosphere 

5.  modeling  of  the  Earth's  dayglow,  Rayleigh  scattering  of  sunlight,  nightglow,  and  auroral 
emissions 

6.  development  of  retrieval  algorithms  to  extract  O2  and  O3  vertical  concentration  profiles  from 
solar  occultation  data  obtained  by  the  SUSIM  instrument  on  UARS 

7.  development  of  software  for  generating  global  images  of  thermospheric  emissions 

8.  modeling  of  miUimeter-wave  emissions  of  the  middle  atmosphere  and  development  of 
correspondirg  retrieval  algorithms 

9.  development  of  a  trace  constituent  climatology  of  the  middle  atmosphere 

10.  improvements  to  (TPFs  optical  emissions  models. 

CPIs  research  efforts  were  carried  out  by  physicists  on-site  at  NRL  and  in-house  at  CPI's 

headquarters  in  Annandale,  VA.  Table  1  identifies  the  on-site  individuals  and  lists  any  key 

experiments  or  models  associated  with  these  researchers  under  this  contract.  Similar  information 

is  provided  in  Table  2  for  in-house  work. 


Table  1.  On-site  CPI  personnel  who  participated  in  this  program.  The  last  two  column  are 
intmded  to  provide  only  an  indication  of  involvement  in  experimental  and  modeling  efforts  by  the 
listed  personnel  over  the  three  year  life  of  this  program. 


Key  Models  Used, 
Developed,  or  Modified 
by  Investieator 


Investigator 

Research  Area 

Satellites  or 
Experiments  Supported 
by  Investigator 

Ms.  K.  A.  Abbott 

Solar  XUV  spectroscopy 

Solar  XUV  images  - 
Skylab 

Dr.  S.  L.  Berg 

Optical  backgrounds 

MSX,  VUE,  IBSS 

Mr.  J.  S.  Morrill 

Dayglow 

RAIDS.  UVLIM. 
MAHRSI 

Mr.  W.  J.  Sawchuck 

Modeling  chemistry  and 
dynamics  of  the  middle 
atmosphere 

MAHRSI 

Dr.  A.  C.  Sterling 

Solar  physics 

Solar-A 

Mr.  S.  E.  Thonnard 

Dayglow,  nightglow 

RAIDS.  SSULI. 
HIRAAS 

PEGFAC 


KINEARTH, 
2D  Chemistry  and 
amics  Model 


Key  to  selected  abbreviations 


MSX 

VUE 

IBSS 

RAIDS 

SSULI 

UVLIM 

MAHRSI 

fflRAAS 

Mf2dsls 

PEGFAC 

KINEARTH 


Midcourse  Space  Experiment 

Visible  Ultraviolet  Experiment 

Infrared  Background  Signature  Survey 

Remote  Atmospheric  and  Itmospheric  Detection  System 

Special  Sensor  Ultraviolet  Limb  Ima^r 

Ultraviolet  Limb  Imager 

Middle  Atmosphere  High  Resolution  Spectrographic  Investigation 
High  Resolution  Airglow  and  Auroral  Spectrograph 


CPI/NRL  dayglow  model  (photoclectron  g-factorsl 
NRL  mesosphere/stratosphere  ID  chemical/dynamics  model 


2 


Table  2.  In-house  CPI  personnel  who  participated  in  this  program.  The  last  two  column  are 
intentted  to  provide  only  an  indication  of  involvement  in  experimental  and  modeling  efforts  by  the 
listed  personnel  over  the  three  year  life  of  this  program 


Investigator 

Research  Area 

Experiments  Supported 
by  Investigator 

Key  Models  Used, 
Developed,  or  Modified 
by  Investigator 

Dr.  D.  E.  Anderson 

Aeronomy 

UVPI 

LOWTRAN,  RAYLIMB 

Mr.  R.  P.  Barnes 

Atmospheric  optical 
emisaons 

MAHRSI 

RAYLIMB 

Dr.  C.-S.  Chang 

Atomic  physics 

SUSIM,  POAM 

Mr.  R.  J.  Cox 

Atmospheric  optical 
emissions 

DE-1  (UV  images), 
UVUM 

Mr.  J.  S.  Evans 

Atmospheric  optical 
emissions 

UVPI,  SSULI 

PEGFAC,  DSDA,  LIMB 

Ms.  D.  L.  Kriebel 

Remote  sensing  of 
middle  atmosphere 

MAS, 

WVMS 

META 

Dr.  J.  D.  Lumpe 

Inversion 

techniques 

SUSIM,  POAM, 
MAHRSI 

LOWTRAN,  SOLOC, 
other  inversion  models 

Dr.  D.  J.  Strickland 

Aeronomy 

UVPI,  SSULI,  VUE 

PEGFAC,  DSDA,  LIMB 

Key  to  selected  abbreviations  not  included  in  Table  1: 


Experiments 

DE-1 

MAS 

POAM 

SUSIM 

UVPI 


Oynamics  Explorer- 1 
Millimeter-wave  Atmospheric  Sounder 
Polar  Ozone  and  Aerosol  Measurement 
Solar  Ultraviolet  Spectral  Irradiance  Monitor 
Ultraviolet  Plume  Listtument 


Models 

DSDA 

LIMB 

RAYLIMB 

LOWTRAN 

SOLOC 

META 


CPl/NRL  photon  transport  code  used  to  calculate  source  functions  of  optically 
thick  lines  in  the  thermosphere  and  ionosphere 

CPI  code  used  to  calculate  limb  radiances  of  dayglow  and  nightglow  emission  features 
CPI  code  which  calculates  limb  profiles  of  Rayleigh  scattered  sunlight 
AFGL  transmission/Rayleigh  scattering  code 

CPI  code  used  to  invert  extinction  proHles  of  full  solar  disk  irradiance  observed 
by  the  NRL  experiment  SUSIM  onboard  satellite  UARS 
CPI/NRL  code  which  calculates  microwave  emissions  of  nrinor  species 
densities  in  middle  and  lo>ver  atmosphere 
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In  summary,  this  contract  provided  the  opportunity  to  conduct  research  on  several  topics 
in  the  general  areas  of  upper  atmospheric  physics,  middle  atmospheric  physics,  and  solar  physics. 
It  also  provided  the  opportunity  to  support  future  experimental  programs  within  NRL's  Space 
Physics  Division  such  as  RAIDS,  SSULI,  SUSIM,  MAHRSI,  and  POAM.  Support  was  in  the 
form  of  laboratory  work,  data  simulations,  and  development  of  algorithms  for  retrieving 
geophysical  information  from  the  listed  instruments.  On  the  SSULI  task,  e.g.,  extensive 
calculations  of  the  dayglow  from  834  A  to  1800  A  were  made  to  estimate  signal  strengths  and 
investigate  the  sensitivity  of  O*  834  A  to  O*  in  the  F-region.  Under  the  SUSIM  task,  an  entirely 
new  retrieval  algorithm  was  developed  which  performs  an  inversion  of  the  occultation  data  to 
(tetermine  number  density  profiles  in  a  model  independent  manner.  In  addition,  a  great  deal  of 
work  went  into  the  development  of  an  accurate  forward  model  code  which  has  been  used 
extensively  in  the  pre-flight  phase  to  generate  realistic  simulations  of  the  SUSIM  occultation  data 
sets.  Under  the  task  involving  remote  sensing  at  millimeter  wavelengths,  the  META  code  has 
developed  to  the  point  where  it  provides  accurate  simulations  of  the  ground  based  millimeter 
wave  emission  measurements  being  performed  by  NRL.  This  is  a  valuable  resource  in  terms  of 
experiment  planning  and  the  development  of  retrieval  schemes  for  this  program.  More  detailed 
discussions  of  these  tasks  and  others  follow  in  the  appendices. 

Papers  have  been  written  which  document  some  of  the  research  performed  in  this  contract. 
A  list  of  these  papers  follows.  The  list  also  includes  abstracts  of  presentations  made  at  AGU 
meetings.  Some  of  them  also  appear  in  the  appendices. 

Lumpe,  J.  D.,  C.  S.  Chang,  and  D.  J.  Strickland,  Atmospheric  Constituent  Density  Profiles  from 
Full  Disk  Solar  Occultation  Experiments,  J.  Quant.  Spectrosc.  Radiat.  Transfer,  46, 483,  1991. 

Strickland,  D.  J.,  R.  J.  Cox,  R.  P.  Barnes,  L.  J.  Paxton,  R.  R.  Meier,  and  S.  E.  Thonnard,  A 
Model  for  Generating  Global  Images  of  Emission  from  the  Thermosphere,  to  be  submitted  to 
Applied  Optics. 

Evans,  J.  S.,  D.  J.  Strickland,  D.  E.  Anderson,  Jr.,  R.  R.  Conway,  H.  W.  Smathers,  D.  F.  Bakeris, 
R.  A.  Swanson,  T.  A.  Chubb,  and  J.  G.  Cardon,  UVPI  Imaging  from  the  LACE  Satellite;  Analysis 
of  Broad  Band  UV  Auroral  Limb  Data,  to  be  submitted  to  J.  Geophys.  Res. 

Morrill,  J.,  W.  Benesch,  Plasma  Preconditioning  and  the  Role  of  Elevated  Vibrational 
Temperatures  in  Production  of  Nj  Excited  State  Vibrational  Distributions,  J.  Geophys.  Res.,  95, 
7711, 1990. 

Carragher,  B.  A.,  J.  Morrill,  W.  M.  Benesch,  Gas  Phase  Molecular  Energy  Transfer  Studies  via 
Time-Resolved  Spectroscopy,  y.  Optic.  Soc.  Am.,  8,  123, 1991. 

Morrill,  J.,  W.  Benesch,  and  K.  Widing,  Electron  Temperatures  in  a  Pulsed  Electric  Discharge  in 
Nj  and  the  Associated  Electron  Excitation  Rate  Coefficients,  7.  Chem.  Phys.,  94, 262, 1991. 

Morrill,  J.,  R.  P.  McCoy,  and  C.  Meng,  Diurnal  Variation  of  Thermospheric  Odd  Nitrogen 
Observed  from  Low  Earth  Orbit,  EOS,  Trans.  Amer.  Geophys.  Union,  70(44),  1236,  1989. 


Morrill,  J.,  D.  Prinz,  and  R.  McCoy,  The  Performance  of  the  MUV  and  NUV  Instruments  in  the 
RAIDS  Experiment,  EOS,  Trans.  Amer.  Geophys.  Union.  71(44),  1486,  1990. 

Fritz,  G.  G.,  D.  K.  Piinz,  R.  P.  McCoy,  R.  R.  Meier,  and  J.  S.  Morrill,  A  Far-Ultraviolet  Imaging 
Spectrograph  for  RAIDS,  EOS,  Trans.  Amer.  Geophys.  Union,  71(44),  1486  1990. 

Newbum,  D.  R.  Goodel,  J.  Morrill,  and  W.  Benesch,  Afterglow  Emission  of  Auroral  Molecular 
Nitrogen  Bands,  EOS,  Trans.  Amer.  Geophys.  Union,  71(16),  213, 1991. 
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ABSTRACT 


SouthOTi  auroral  limb  images  have  bew  recorded  by  the  wide  angle  UV  camera  within 
the  Ultraviolet  Plume  Instrument  (UVPI)  onboard  the  Low-power  Atmospheric  Compensation 
Experiment  (LACE)  satellite.  Tangent  altitudes  within  each  image  have  been  determined  using 
recorded  stars  along  with  the  known  position  of  the  satellite.  Range  to  the  observed  aurora  has 
been  determined  by  stereoscopic  techniques  using  images  of  die  same  auroral  forms  collected 
over  several  seconds.  An  auroral  model  has  been  used  to  fit  two  limb  profiles  taken  from  the 
same  image.  The  parameters  varied  within  the  model  to  achieve  the  given  fits  are  1)  number  of 
auroral  forms  within  the  observing  region,  2)  their  locations,  3)  thicknesses,  and  4)  spectral 
characteristics  of  the  precipitating  electron  spectra.  The  model  provides  altitude  profiles  of 
volume  emission  for  the  relevant  features  based  on  electron  transport  calculations  starting  with 
a  user  specified  precipitating  electron  spectrum.  Volunw  emission  for  multiple  precipitation 
patterns  is  then  distributed  within  the  observing  region  followed  by  line-of-sight  integrations 
producing  limb  profiles.  Excellent  fits  to  the  entire  observed  profiles  were  obtained  using  a 
combination  of  diffuse  aurora  and  a  single  discrete  aurora  placed  near  the  tangent  point.  The 
fitting  process  suggests  that  a  good  quantitative  description  of  the  aurora  has  been  achieved  in 
terms  of  precipitating  electron  spectral  characteristics  and  spatial  distribution  of  this  precipitation. 


1.  INTRODUCTION 


The  subject  of  this  paper  is  interpretation  of  a  sequence  of  UV  images  acquired  during 
a  limb  scan  of  a  southern  aurora  by  the  UVPI  (T^traViolet  Plume  jhstrument)  experiment  on 
satellite  LACE  (Low-power  Atmospheric  Compensation  Experiment).  Spatial  and  energy 
characteristics  of  electron  precipitation  in  the  observed  aurora  are  inferred  from  a  detailed  fitting 
of  calculated  broad  band  UV  limb  profiles  to  the  data.  To  our  knowledge,  this  is  the  first  paper 
to  address  satellite  observed  aurora  on  the  limb  in  a  quantitative  fashion  without  the  use  of 
tomographic  inversion  or  algebraic  reconstruction  techniques  (see  e.g.,  Solomon  et  al.,  1984, 
Solomon  et  al.,  1985,  McDade  and  Llewellyn,  1991,  and  references  therein).  Past  satellite  UV 
imaging  experiments  have  obtained  auroral  images  containing  mostly  disk  emission.  Any  auroral 
limb  data  currently  available  have  either  been  on  much  coarser  tangent  altitude  scales  or  were 
measured  tomographically  (i.e.  limb  scans  through  one  or  more  emitting  regions  observed  from 
different  locations).  Satellite  observations  include  those  from  Dynamics  Explorer  [Frank  and 
Craven,  1988],  Viking  (Anger  et  al.,  1987],  HILAT  [Meng  and  Huffman,  1984],  and  PolarBEAR 
[Schenkel  et  al.,  1986]  (also  see  Solomon  et  al.,  1988,  Solomon,  1989,  McDade  et  al.,  1991, 
Swenson  et  al.,  1990). 

Our  interest  is  to  understand  the  spectral  content  of  the  auroral  emission  being  observed 
and  to  investigate  how  well  spatial  and  energy  characteristics  of  the  aurora  can  be  inferred  from 
single  scan  limb  observations.  Four  characteristics  of  the  aurora  determine  its  variation  in 
luminosity  with  tangent  altitude.  They  are  1)  number  of  forms  along  the  line-of-sight,  2)  their 
locations,  3)  their  horizontal  extents,  and  4)  hardnesses  and  energy  content  of  the  incident 
electron  spectra.  We  will  quantitatively  address  the  relationship  between  luminosity  variations 
across  the  limb  and  these  properties  using  the  auroral  model  of  Strickland  et  al.  [1976]  (also  see 
Strickland  et  al.,  1983,  Strickland  et  al.,  1989,  Basu  et  al.,  1993,  and  Strickland  et  al.,  1993). 
The  model  is  then  used  to  fit  two  limb  profiles  taken  from  a  single  tracker  camera  image.  As 
will  be  shown,  this  requires  a  minimum  of  two  auroral  forms  along  the  line-of-sight.  While  the 
above  inferred  properties  are  not  unique,  we  believe  that  the  four  parameters  identified  above  can 
be  sufficiently  restricted  to  provide  quantitative  information  on  spatial  and  energy  characteristics 
of  the  aurora.  — 


-  1  - 


2.  THE  LACE  SATELLITE 


The  LACE  satellite  was  designed  and  built  by  the  Naval  Research  Laboratory  (NRL)  for 
the  Strategic  Defense  Initiative  Organization.  TTie  satellite’s  primary  mission  was  to  provide  a 
sensing  array  in  orbit  for  evaluating  the  effect  of  active  compensation  of  a  ground  based  laser 
system  on  light  propagated  through  the  atmosphere.  On  February  14,  1990,  the  LACE  satellite 
was  launched  into  a  circular  orbit  of  540  km  with  an  inclination  of  43“.  The  LACE  spacecraft 
was  designed  to  operate  the  primary  experiment  for  a  30  month  lifetime.  It  is  a  gravity-gradient 
stabilized  satellite  with  attitude  measurement  curacy  better  than  1“.  To  support  the  LACE 
spacecraft,  there  were  three  ground  stations  operated  by  NRL,  one  permanent  and  two 
transportable.  LACE  was  decommissioned  on  February  14,  1993. 

The  UVPI  was  a  pointable,  telescopic  instrument  mounted  on  the  Earth  facing  side  of  the 
LACE  spacecraft  Its  mission  was  to  collect  images  of  rocket  plumes  in  the  near-  and  mid¬ 
ultraviolet  regimes  from  a  space  based  platform.  It  was  also  used  to  collect  background  image 
data  of  the  Earth,  Earth’s  limb,  and  celestial  objects.  Background  imaging  performed  by  UVPI 
included  the  day,  night,  and  dawn  limb,  aurora,  sunlit  and  moonlit  clouds,  and  the  Earth’s 
hardbody. 

The  UVPI  consisted  of  two  intensified  CCD  cameras  that  were  boresighted  and  shared 
a  Maksutov-Cassegrain  telescope.  The  tracker  camera  was  sensitive  in  the  near-ultraviolet^visible 
wavelength  region,  had  a  15  times  larger  field-of-view  than  the  plume  camera,  and  was  used  to 
locate,  acquire,  and  track  a  target  The  plume  camera  collected  images  at  near-  and  mid- 
ultraviolet  wavelengths  selectable  with  a  filter  wheel.  Fig.  1  shows  a  schematic  of  the  intensifred 
camera.  Table  1  lists  characteristics  of  the  optics  and  each  camera.  The  camera  images  were 
digitized  and  transmitted  to  the  ground  or  were  recorded  onboard  for  later  transmission. 

The  optics  for  both  camera  systems  are  shown  in  Fig.  2.  The  UVPI  was  mounted  within 
the  satellite  and  looked  through  an  aperture  in  the  Earth  oriented  end  panel.  By  use  of  a 
gimballed  mirror  the  UVPI  had  a  field  of  regard  of  a  50“  half-angle  cone  about  the  satellite’s 
nadir.  When  not  in  use,  a  door  covered  the  UVPI  aperture.  Attached  to  the  inside  of  this  door 
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was  a  flat  mirror  which,  when  the  door  was  opened  part  way,  allowed  the  UVPI  cameras  to  view 
celestial  objects  or  the  Earth’s  limb.  Because  of  the  hinge  arrangement  of  the  calibration  door, 
limb  viewing  was  only  possible  out  the  right  side  of  the  spacecraft,  which  was  due  south  when 
the  spacecraft  was  moving  due  east  at  43°  south  or  north  latitude.  No  viewing  of  the  northern 
aurora  was  possible. 

The  images  analyzed  in  this  paper  were  acquired  at  a  rate  of  5  s\  The  camera  output 
was  an  analog  signal  read  from  a  CCD  chip  as  a  sequence  of  interlaced  image  fields.  The  analog 
signal  was  immediately  digitized  to  produce  a  shifted  (temporally  delayed)  image  field  once 
every  1/60  s.  Pairs  of  interlaced  field  readouts  were  summed  to  provide  a  full  digitized  image 
frame  once  every  1/30  s.  In  normal  mode,  every  sixth  frame  was  sent  to  the  ground  in  real-time, 
or  was  sent  to  a  tape  recorder  and  subsequently  transmitted  to  a  ground  station.  Each  image 
consists  of  an  array  of  251  x  240  8-bit  numbers  that  measures  the  flux  distribution  on  the  CCD 
sensor  plane. 

The  photometric  range  and  sensitivity  of  the  UVPI  cameras  were  optimized  for  nighttime 
operations.  That  is,  the  instrument  was  designed  to  view  and  track  relatively  bright  targets 
against  a  dark  background.  The  UVPI  tracker  camera  could  image  stars  as  dim  as  visual 
magnitude  7.  The  plume  camera  had  a  commandable  filter  wheel  allowing  selection  from  four 
filters  of  varying  passbands.  The  relative  spectral  responsivity  was  determined  from  transmission 
and  reflectivity  measurements  of  the  optical  components  of  the  system  and  measurement  of  the 
photoelectric  efficiency  of  the  image  intensifier  photocathode.  The  absolute  responsivity  of  the 
system  was  determined  from  post-launch  observations  of  stars.  The  average  deviation  in  the  ratio 
of  observed  to  calculated  star  responses  was  about  10-15  percent  in  the  plume  camera,  and  about 
16  percent  in  the  tracker  camera.  Reynolds  et  al.,  [1989]  give  a  more  complete  description  of 
the  UVPI  and  its  requirements. 

On-orbit  calibration  was  accomplished  by  observing  and  actively  tracking  stars  with  a 
known  spectrum.  Calibration  star  measurements  were  made  frequently,  about  four  to  eight  times 
a  month,  mostly  to  support  UVPI’s  primary  mission  and  to  collect  UV  spectral  data  on  stars. 
Smathers  et  al.,  [1989]  gives  a  description  of  pre-launch  calibration  and  performance.  The  data 
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presented  in  this  paper  come  from  the  tracker  camera  over  the  wavelength  interval  from  250  to 
450  nm.  For  these  data,  the  tracker  camera  was  not  centered  on  the  observed  auroral  forms;  thus 
the  plume  camera  did  not  record  them. 

3.  OBSERVING  CONDITIONS 

The  images  of  interest  were  recorded  during  orbit  2570  on  August  3, 1990  while  viewing 
the  southern  aurora  in  the  morning  sector  at  approximately  17:20  GMT.  The  q>acecraft  was  in 
darkness  over  the  entire  viewing  period  which  provided  optimum  efficiency  of  the  CCD  cameras 
and  associated  optics.  The  observed  aurora  was  also  in  darkness.  The  index  was  3  which 
reflec.o  moderate  geomagnetic  activiQr.  The  spacecraft  was  near  the  southernmost  point  of  the 
orbit  (43.1°  latitude,  110°  E  longitude)  off  the  southern  coast  of  Australia.  Viewing  was  to  the 
south  and  for  the  data  of  interest,  the  horizontal  center  line  of  the  image  occurred  at  about  100 
km  tangent  altitude  with  a  corresponding  slant  range  of  2475  km  where  the  spatial  coverage  of 
the  tracker  camera  was  374  km  x  491  km  (vertical  by  horizontal).  Several  auroral  forms  are 
present  in  the  data  collected.  Some  300  firames  were  recorded  by  the  tracker  camera  during  an 
observation  period  of  100  s.  Plate  1  shows  a  composite  of  selected  ovnlapping  firames  from 
which  the  data  to  be  presented  in  Section  6  were  taken. 


0  4.  TANGENT  ALTITUDE  AND  RANGE  DETERMINATIONS 

The  apparent  height  of  the  aurora  above  the  limb  was  determined  from  the  apparent  height 
of  a  background  star  observed  in  the  same  images,  using  the  star  coordinates  and  the  orbital 
#  elements  of  the  spacecraft  (s/c).  Auroral  range  was  determined  from  the  apparent  backward 

motion  of  the  aurora  (v,„).  This  motion  results  primarily  from  the  forward  orbital  velocity  but 
is  also  affected  by  rotation  of  the  spacecraft  (one  rotation  per  orbit)  and  rotation  of  the  Earth. 
Assuming  a  stationary  aurora,  v^,  is  given  by 

^  ''-r  =  -V  +  X  »-U>S  +  ^  R-r*  “r 
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where  is  the  orbital  velocity,  (O^  =  x  is  the  angular  velocity  of  the  spacecraft  in  its 
orbit,  r^os  is  the  line>of-sight  vector  fix)m  the  spacecraft  to  the  aurora,  R,.,  and  R^  are  vectors 
from  Earth  center  to  the  aurora  and  the  s/c  respectively,  is  the  angular  velocity  of  Earth,  and 
u,  is  the  unit  vector  v^  /  |v^  |.  Solving  for  the  range  r,^ 

,  _  V.,  ^  cos  cos  (5) 

sin(a) 

where  a  is  the  depression  angle  of  r^^,  is  the  latitude  of  the  aurora,  and  5  is  the  angle 
between  and  due  East  For  the  data  of  interest  here,  the  contributions  to  apparent  auroral 
velocity  were  v^  =  -7.204  km/s,  lo)^*  x  r,^|  =  0.992  km/s,  and  |f2g  x  R,^  '  u,|  =  0.229  km/s. 
The  range  to  the  aurora  seen  near  a  tangent  altitude  of  1 10  km  was  r^.  =  2475  km. 


5.  MODELLING  OF  AURORAL  EMISSION  IN  A  LIMB  VIEWING  GEOMETRY 

We  use  the  electron  transport  model  of  Strickland  et  al.  [1976]  generalized  to  three  neutral 
density  constituents  for  calculating  the  electron  flux  (|i(z,E,p)  in  c  cm*V‘eV  ‘sf  ‘  starting  with  a 
precipitating  electron  spectrum  <t>,(E,p)  in  the  same  units.  The  independent  variables  are  altitude 
z  in  km,  energy  E  in  eV,  and  cosine  of  the  pitch  angle  p.  Maxwellian  and  Gaussian  distributions 
isotropic  over  the  downward  hemisphere  have  been  used  to  characterize  ()>„  in  this  work.  We  use 
Gaussian  distributions  to  describe  the  precipitating  electron  spectrum  for  the  discrete  aurora  and 
either  Gaussian  or  Maxwellian  distributions  to  describe  electron  spectrum  for  the  diffuse  aurora. 
High  and  low  energy  tails  are  included  as  given  in  Strickland  et  al.  [1993].  For  either 
representation,  the  key  parameters  are  energy  flux  Q  in  ergs  cm'V*  and  characteristic  enwgy  E,. 
For  Maxwellians,  E„  is  half  the  mean  energy  while  for  Gaussians,  it  refers  to  the  energy  of  the 
Gaussian  peak.  The  approximate  range  of  E„  for  diffuse  aurora  is  O.S  to  3  keV.  For  discrete 
aurora,  a  few  keV  to  about  10  keV  is  common. 


From  (|>(z,E,p),  we  calculate  volume  excitation  and  ionization  rates  with  the  appropriate 
cross  sections  for  use  in  chemistry  modeling  and  in  specifying  emission  rates.  Examples  of 
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excitation  and  ionization  rates  for  a  variety  of  may  be  seen  in  Strickland  et  aL  [1989]. 
Figure  3  shows  emission  rates  for  lines  and  band  systems  within  the  wavelengdi  interval  of 
interest  here  (180-450  nm).  A  Gaussian  distribution  as  described  above  with  Q  »  10  ergs  cm'V^ 
and  Eg  s  S  keV  was  used  to  characterize  the  precipitating  electron  spectrum.  The  emission 
features  include  the  atomic  lines  OI  297.2  nm  and  NI  346.6  nm  and  bands  of  die  following 
molecular  systems:  N,  LBH,  Nj  VK,  Nj  2PG,  Nj*  ING,  Oj  HerzbCTg  I  and  H,  and  O2 
Chamberlain.  The  Oj  volume  emission  rate  profile  does  not  include  nightglow.  A  nightglow 
component  will  be  included  when  comparing  with  UVPI  data.  The  emission  at  346.6  nm  arises 
from  dissociative  excitation  of  N,.  Quenching  is  responsible  for  the  slower  altitude  variation  of 
the  Nj  VK  and  OI  297.2  nm  profiles. 


Figure  4  illustrates  the  spectral  characteristics  of  limb  emission  from  N,  LBH,  N,  VK, 
Nj  2PG,  and  Nj*  ING.  The  incident  electron  spectnim  is  the  same  as  used  to  generate  the 
volume  emission  rates  in  Figure  3.  The  arc  thickness  is  100  km  with  a  position  of  *1000  km. 
The  tangent  altitude  is  -25  km. 

Figure  5  shows  a  sketch  illustrating  how  slant  intensities  for  various  auroral  forms  were 
computed.  The  narrow  shaded  region  represents  a  discrete  arc  at  the  tangent  point  while  the  wide 
shaded  region  represents  diffuse  aurora  between  the  tangent  point  and  the  spacecraft  We  assume 
a  given  thickness  for  an  arc  such  as  100  km  and  assume  the  precipitating  electron  q)ectrum  to 
be  the  same  over  this  thickness.  The  arc  is  then  placed  at  some  position  along  the  line-of-sight 
as  shown  in  the  figure.  The  emitting  layer  shown  in  Figiue  5  is  meant  to  represent  the  vertical 
^''tent  of  the  auroral  arcs.  The  emitting  region  generally  begins  around  80  km  and  ends  at  about 
300  kin  with  peak  intensities  occurring  in  the  vicinity  of  1 10  km.  The  arc  position  is  referenced 
to  the  tangent  point.  For  example,  a  distance  of  -5(X)  km  would  refer  to  an  arc  whose  center  is 
5(X)  km  from  the  tangent  point  in  the  direction  of  the  observer. 


Figure  6  provides  examples  of  limb  intensity  profiles  of  the  ING  system  in  Rayleighs 
for  a  100  km  thick  arc  placed  at  the  following  locations:  -10(X),  -500,  0,  and  500  km.  The 
incident  electron  spectrum  is  the  same  one  used  to  generate  the  results  in  Hgures  3  and  4.  As 
expected,  the  profile  for  the  arc  at  the  tangent  point  peaks  at  a  higher  tangent  altitude  than  those 
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in  the  near-  and  far-fiekl.  Its  maximum  value  occurs  at  1  IS  km  which  is  close  to  the  altitude  of 
maximtun  volume  emission.  The  4-1000  km  profile  has  the  sharpest  fall  off  on  its  bottom  side 
due  to  geometry  considerations.  Attenuation  along  the  line-of-sight  was  not  taken  into  account 
in  the  calculations.  It  is  not  important  for  the  given  system  except  possibly  at  tangoit  altitudes 
well  below  100  km  for  arcs  on  the  far  side  of  the  tangent  point  In  this  case,  Rayleigh  scattering 
will  reduce  the  intensity  on  the  bottom  side  of  the  profiles  but  not  enough  to  significantly  change 
their  shr^.  For  emission  below  300  nm  (see  Figure  4),  an  additional  source  of  attenuation  is 
ozone.  Again,  it  is  not  important  except  for  emission  originating  on  the  far  side  of  the  tangent 
point  and  being  recorded  at  tangent  altitudes  well  below  100  km.  For  the  data  to  follow,  we  will 
argue  that  most  of  the  recorded  emission  at  small  tangent  altitudes  comes  from  the  near  side  of 
the  tangent  point  and  thus  experiences  little  or  no  attenuation  by  either  absorption  or  scattering. 

Figure  7  shows  the  sensitivity  of  limb  intensities  to  using  Gaussian  distributions  with 
tails  and  for  a  Q  value  of  10  ergs  cm'^  s'*.  The  arc  has  a  thickness  of  100  km  with  its  center 
located  at  the  tangent  point.  As  in  Figure  6,  we  have  selected  the  ING  system  for  displaying 
the  results.  The  same  variation  versus  E,  as  well  as  versus  arc  location  (previous  figure)  occurs 
for  any  other  emission  feature  within  the  180  to  450  nm  region.  A  comparison  of  Figure  5  and 
Figure  6  shows  that  changes  in  arc  location  mimic  changes  in  E^.  As  an  arc  moves  away  fix>m 
the  tangent  point,  the  limb  intensity  profile  peak  decreases  in  tangent  altitude  as  if  an  increase 
in  E,  had  occurred.  Unlike  the  latter  case,  however,  the  tangent  altitude  of  the  peak  can  decrease 
to  values  considerably  below  100  km  which,  if  interpreted  in  terms  of  E„  would  lead  to 
unrealistically  large  values  in  this  parameter.  It  is  important  to  be  able  to  distinguish  between 
changes  in  the  limb  intensity  profile  due  to  changes  in  arc  position  versus  changes  in  E^. 
Otherwise,  unique  information  on  the  spatial  distribudon  of  emission  and  characteristics  of  the 
associated  precipitation  cannot  be  obt?in^d  from  limb  intensity  profiles.  In  the  next  section,  we 
will  argue  that  it  is  possible  to  distinguish  between  these  two  effects  using  limb  data  that  provide 
a  well  defined  shape  to  the  limb  intensity  profile.  Support  for  tliis  argument  comes  from  Figures 
6  and  7  where  it  is  seen  that  considerably  more  change  in  shape  (ignoring  where  the  tangent 
altitude  of  the  peak  occurs)  occurs  as  changes  compared  to  changes  in  arc  location. 


6.  FITTING  MODEL  LIMB  PROFILES  TO  THE  DATA 


The  shape  of  an  auroral  limb  profile  is  a  function  of  the  number  of  aurtmd  forms  along 
the  line-of-sight  and  individual  characteristics  of  these  foims.  These  characteristics  fall  into 
spatial  and  spectral  categories.  Within  the  spatial  category  are  thickness  and  location.  Within 
the  spectral  category  are  Q,  E,,  and  the  type  of  distribution  (Maxwellian  or  Gaussian  in  this 
work).  This  totals  to  six  parameters/descriptions  which  raises  the  question  of  uniqueness  when 
a  good  fit  to  the  data  is  achieved.  We  will  not  address  this  issue  in  any  detail  here  but  will  show 
that  real  constraints  can  be  placed  on  their  infened  values  through  observing  the  quality  of  fit 
as  these  parameters/descriptions  are  changed.  An  example  of  one  condition  imposed  by  the  given 
data  is  that  at  least  two  auroral  forms  are  required  to  fit  the  data.  A  good  fit  could  not  be 
achieved  using  one  auroral  structure  regardless  of  its  spatial  and  spectral  characteristics.  More 
will  be  said  about  this  below.  Limb  profile  shape  is  most  sensitive  to  the  location  of  the  auroral 
form  and  to  E^.  The  degree  of  sensitivity  was  illustrated  in  Hgures  6  and  7  as  discussed  in  the 
previous  section.  Sensitivity  to  the  remaining  parameters/descriptions  is  less  important  but  all 
parameters  were  varied  during  the  data  fits  to  be  discussed  next 

Figure  8  shows  the  best  fits  to  two  data  sets  taken  from  the  same  tracker  camera  image. 
As  noted  in  Section  2,  the  data  were  collected  on  pass  2570  looking  south  at  the  southern 
nighttime  aurora  on  August  3,  1990  (see  plate  1).  The  error  bars  shown  with  the  data  represent 
±  1  standard  deviation  (assuming  Poisson  statistics).  We  will  refer  to  the  data  in  the  upper 
Oower)  panel  as  Set  #1  (Set  #2)  data.  The  two  measured  profiles  possess  similar  peaks  and 
topside  scale  heights.  The  bottomsides,  however,  are  different  and  will  be  further  discussed 
below.  The  model  profiles  each  comprise  a  nightglow  (Oj  Herzberg  I,  n  and  Chamberlain  band 
systems)  and  two  auroral  components.  The  nightglow  component  used  in  this  analysis  is  based 
on  the  work  of  Stegman  and  Murtagh  [1991]  and  is  consistent  with  the  results  obtained  therein. 
For  the  brightness  of  the  given  aurora,  the  nightglow  component  is  unimportant. 

As  mentioned  above,  a  good  fit  could  not  be  achieved  with  a  single  auroral  form. 
Excellent  fits  were  achievable,  however,  by  introducing  two  auroral  forms,  one  near  the  tangent 
point  and  the  other  in  front  of  it.  Table  2  provides  details  by  identifying  the  six 
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parameters/conditions  discussed  above.  The  inferred  location  of  the  auroral  forms  provided  in 
Table  2  refer  to  distance  from  the  tangent  point  towards  the  observer.  The  incliMled  error 
estimates  are  accqrtable  deviations  from  the  infened  values  and  reflect  the  sensitivity  analysis 
performed  on  the  data.  For  values  outside  these  deviations,  the  tits  to  the  data  become  degraded. 
It  should  be  noted  that  any  errors  in  calibration  are  directly  reflected  in  the  magnitudes  of  the  Q 
values. 


For  either  data  set,  an  arc  at  the  tangent  point  is  needed  to  tit  the  peak  region.  The 
inferred  thickness  is  300  km  and  suggests  a  discrete  aurora.  A  second  auroral  structure  in  front 
of  the  tangent  point  is  required  to  tit  the  topside  and  bottomside.  If  placed  behind  the  tangent 
point,  a  tit  to  the  bottomside  cannot  be  achieved.  The  inferred  thickness  is  approximately  1000 
km  for  either  data  set  and  suggests  diffuse  aurora.  The  inferred  £„’$  for  the  second  auroral 
structure  are  typical  hardnesses  of  a  diffuse  aurora.  The  difference  in  bottomsides  of  the  two  data 
sets  is  due  to  differences  in  the  diffuse  aurora  in  front  of  the  tangent  point  To  achieve  the 
greater  falloff  in  Set  2.  a  somewhat  narrower  diffuse  region  closer  to  the  observation  point  was 
needeo  as  quantified  in  Table  2. 


7.  DISCUSSION  AND  CONCLUSIONS 

A  detailed  analysis  of  single  scan  satellite  observed  auroral  limb  emission  in  terms  of 
structure  and  spectral  characteristics  has  not  been  performed  before.  A  general  impression  witiiin 
the  tield  has  been  that  little  useful  quantitative  information  comes  from  such  data  given  the 
general  lack  of  knowledge  about  structure  along  the  line-of-sighL  This  analysis  suggests  that  an 
accurate  measurement  of  limb  emission  from,  say,  0  to  200  km  in  tangent  altitude  allows  one  to 
reasonably  specify  the  gross  structural  characteristics  of  the  emission.  For  example,  it  seems 
clear  from  this  analysis  that  two  structured  forms  with  the  characteristics  given  in  Table  2 
produced  the  observed  limb  profiles.  It  then  follows  that  one  can  also  infer  spectral 
characteristics  of  the  precipitating  electrons  as  was  done  here.  The  shortcoming  of  the  present 
results  is  that  we  did  not  apply  an  error  propagation  model  to  produce  good  quantitative  error 
bars  which  would  address  the  issue  of  uniqueness.  A  first  step  was  made  -in  the  form  of  a 


sensitivity  study  in  which  the  quality  of  fit  was  monitored  as  parameter  values  were  varied.  The 
next  step  is  to  perform  similar  analyses  of  limb  data  for  other  auroras  exhibiting  distinctly 
differmt  spatial  and/or  spectral  characterisdcs.  This  will  jnovide  more  information  on  constraints 
imposed  by  the  fitting  process.  Independent  information  on  spatial  structure  fitmi  ground  based 
or  space  based  unagery  would  also  be  useful. 

In  conclusion,  we  believe  that  Table  2  provides  a  good  overall  description  of  the  aurora 
producing  the  data  in  Sets  1  and  2  based  on  our  detailed  examination  of  fitting  errras  as  the 
parameters  listed  in  Table  2  were  varied. 
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PLATES 


Plate  1.  Tracker  camera  conqwsite  limb  image  of  the  southern  aurora  taken  from  -17:20  to 
17:22  GMT  on  August  3, 1990.  The  horizontal  center  line  of  the  image  corresponds  to  a  tangent 
altitude  of  -100  km.  The  white  dots  show  the  centers  of  the  individual  frames  used  to  make  the 
composite.  Where  two  frames  overlap,  data  are  taken  from  the  most  recently  recorded  frame. 
The  frame  centers  are  connected  to  emphasize  the  motion  of  the  canneras  (tracker  and  plume). 
The  vertical  motion  is  produced  by  a  moving  mirror  while  that  in  the  horizontal  direction  is 
produced  by  the  motion  of  the  satellite.  The  data  to  be  analyzed  were  taken  firom  two  vertical 
cuts  through  the  image.  Data  set  1  (set  2)  contains  limb  emission  from  the  Inightest  region  of 
the  structure  on  the  left  (right). 


HGURES 

Figure  1.  Schematic  diagram  of  the  CCD  cameras  of  the  Ultraviolet  Plume  Instrument 

Figure  2.  Diagram  of  the  shared  optics  of  UVPI  including  Maksutov-Cassegrain  telescope. 

Figure  3.  Volume  emission  rates  of  the  key  features  contributing  to  the  spectral  region  from  180 
to  450  nm.  The  applied  incident  electron  spectram  is  a  Gaussian  distribution  with  low  and  high 
energy  tails  having  Q  and  E,  values  of  10  ergs  cm'^  s'*  and  5  keV,  respectively.  The  profile 
labeled  "HBRG"  gives  emission  from  the  O,  Herzberg  I,  Herzberg  n,  and  Chamberlain  systems 
by  electron  impact 

Figure  4.  Spectra  for  individual  band  systems  in  Rayleighs/nm  seen  at  a  tangent  altitude  of  -25 
km  from  180  to  450  nm.  The  arc  used  to  generate  these  results  is  at  -1000  km  and  has  a 
thickness  of  100  km.  The  O,  Herzberg  I,  Herzberg  II,  and  Chamberlain  band  systems  are  not 
displayed  for  lack  of  synthetic  spectral  models  at  this  time.  Their  total  emission,  however,  is 
included  in  the  figure. 
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Figure  S.  Schematic  diagram  for  illustrating  how  auroral  structure  is  treated  in  the  calculations. 

Figure  6.  Limb  intensity  profiles  of  N/  ING  emission  for  a  100  km  thick  arc  at  four  locations 
along  the  line-of-sight  The  locations  are  relative  to  the  tangent  point  The  incident  electron 
spectrum  is  the  same  used  to  goierate  the  results  in  Figure  3.  The  observation  altitude  is  SOO 
km. 


Figure  7.  Limb  intensity  profiles  of  ING  emission  versus  for  Gaussian  distributions  with 
tails.  The  energy  flux  Q  =  10  ergs  cm*V^  The  arc  is  centered  at  the  tangent  point  with  a 
thickness  of  100  km.  The  observation  altitude  is  SOO  km. 

Figure  8.  Fits  to  two  auroral  limb  profiles  taken  from  a  single  UVPI  tracker  camera  image. 
Each  of  the  two  calculated  profiles  is  composed  of  the  three  components  displayed  in  the  figtue. 
For  the  given  brightness  of  this  aurora,  the  nightglow  component  is  unimportant  The 
observation  altitude  is  535  km. 
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Tables 


Table  1.  UVPI  tracker  and  plume  camera  characteristics  and  performance. 


Parameter 

Plume  Camera 

Tracker  Camera 

Spectral  region 

195-350  nm 

255-450  nm 

Number  of  filters 

4 

1 

Time  for  filter  change 

1.7  sec 

n/a 

Frame  rate,  normal 

5  per  sec 

5  per  sec 

Frame  rate,  zoom* 

30  per  sec 

30  per  sec 

Super  pixels,  normal 

251  X  240 

251  X  240 

Super  pixels,  zoom* 

91  X  112 

91  X  112 

Digitization 

8  bits/pixel 

8  bits^ixel 

Image  data  rate 

2.5  MB/sec 

2.5  MB/sec 

Reid  of  view 

0.184  X  0.137  deg 

2.60  X  1.98  deg 

Field  of  regard 

100  X  97  deg 

100  X  97  deg 

Pixel  field  of  view 

12.8  X  10.0  pR 

181  X  143.9  pR 

System  resolution 

80  -  100  pR 

220  -  250  pR 

Pixel  exposure  time 

33  ms 

^3  ms  (adjustable) 

Frames  integrated 

1  -6 

n/a 

Telescope  aperture 

10  cm 

10  cm 

♦reduced  field  of  view 
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Table  2.  Parameters  used  to  fit  the  data  in  Bgure  8.  M/G  in  the  third  column  r^ers  to  Maxwdlian 
or  Gaussian  distributions. 


Data 

Set 

Auroral 

Structure 

(M^) 

m 

Q 

(ergs  cm  *  s  ‘) 

Thickness 

(km) 

Location 

Oon) 

#1 

#1 

G 

4.5  ±  .5 

8.2  ±  2.5 

300  ±  100 

0150 

#1 

#2 

G 

2.2  ±  .3 

3.8  ±  1.3 

1  100  ±250 

-550  1  50 

#2 

#1 

G 

7.0  ±  1 

4.2  ±  1.4 

300  ±  100 

0150 

#2 

#2 

M 

1.0  ±  .5 

5.7  ±  1.8 

9501250 

-350  1  50 
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1.  INTRODUCTION 


# 


# 


In  this  report,  we  begin  to  assess  how  well  the  daytime  F-region  O^  density  profile 
(effectively,  the  F-region  electron  density  profile)  can  be  inferred  firom  O^  834  A  limb 
profiles  to  be  obtained  by  the  EUV/FUV  limb  sensor  SSULI  (Special  Sensor  Ultraviolet  I^b 
Imager)  on  future  DMSP  (^fense  Meteorogical  Satellite  J^gram)  SD3  satellites.  Emphasis 
will  be  on  the  shape  of  the  834  A  limb  profile  and  how  sensitive  it  is  to  changes  in  the  O^ 
concentration  which  we  will  characterize  by  Chapman  functions.  The  independent  variables 
are  peak  density  Ng^x,  altitude  of  the  peak,  h,,,,,  and  topside  scale  height  Hiop.  Changes  in  the 
834  A  limb  profile  will  be  examined  for  changes  in  N„,„  firom  5x10*  to  6x10*  cm'®  and  hg„j 
from  300  to  600  km.  The  topside  scale  height  is  kept  fixed  at  120  km.  The  dependence  of 
834  A  emission  on  O*  is  through  the  multiple  scattering  of  its  photons  which  are  initially 
produced  by  photoionization  of  O  with  minor  contributions  from  electron  impact  on  O  and 
resonance  scattering  of  834  A  solar  photons.  The  dayglow  model  PEGFAC  (Strickland  and 
Meier,  1982)  is  used  to  calculate  the  initial  834  A  volume  emission  rates  followed  by  use  of 
code  DSDA  (Strickland  and  Anderson.  1983)  to  calculate  the  corresponding  final  rates 
obtained  from  multiple  scattering  calculations.  Code  LIMB  is  then  used  to  perform  line-of- 
sight  integrations  giving  834  A  limb  intensity  profiles  in  units  of  Rayleighs.  For  this  study, 
the  solar  zenith  angle  (SZA)  was  held  fixed  along  the  line-of-sight  We  refer  the  reader  to 
the  studies  of  Anderson  and  Meier  (1985),  McCoy  et  al.  (1985),  Cleary  et  aL  (1989),  Paxton 
and  Strickland  (1988),  Strickland  et  al.  (1993),  and  Link  et  al.  (1993)  for  previous  woric  on 
this  problem. 


2.  KEY  INPUT  PARAMETERS 


The  model  atmosphere  used  to  calculate  the  834  A  dayglow  is  MSIS-86  (Hedin, 

1987).  Table  1  shows  the  assumed  values  of  its  input  parameters.  We  have  chosen  an  Fj,,, 
value  of  200  which  corresponds  to  moderate  to  high  solar  activity.  This  value  is  also  used  to 
specify  the  solar  EUV  spectrum  which  wilt  be  discussed  below.  The  resonance  and  pure 
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Table  1. 


Values  of  MSIS  input  parameters  used  to  calculate  834  A  results. 


Model  Atmosphere 


Model 

MSIS-86 

Year 

1990 

Day 

81 

UT(sec) 

0.0 

Geo.  lat 

67.6“ 

Geo.  Ion 

224.1“ 

^10.7 

200 

<Pl0.7> 

200 

Ap 

40 

Table  2.  Pure  absorption  and  resonance  cross  sections  of  the  834  A  triplet  lines  used  to 
perform  the  multiple  scattering  calculations. 


Waveleneth  (A)  Pure  Absorption  Cross  Section  (cm^ )  Resonance  Cross  Section  (cm^l 


N, 

0* 

O 

832 

1.00x10*’ 

3.30x1(7*’ 

4.20x1(7** 

4.14x1(7’ 

833 

4.00x10*’ 

1.40x1(7*’ 

4.20x1(7** 

8.80x1(7’ 

834 

1.20x1(7** 

l.lOxlO-*’ 

4.20x1(7** 

1.22x1(7’ 

-2- 


absorption  cross  sections  needed  to  perform  the  multiple  scattering  calculations  are  shown  in 
Table  2.  The  solar  EUV  spectrum  comes  from  Hinteregger  (1981a  and  1981b)  and  is  shown 
in  Figure  1.  Hinteregger  provides  an  empirical  model  for  wavelength  dependent  scaling  of 
his  reference  spectrum  to  account  for  solar  activity.  The  parameter  used  for  the  scaling  is 
^10.7-  As  noted  above,  we  have  used  a  value  of  200  for  producing  this  spectrum. 


3.  SAMPLE  DMSP  ORBIT  AND  MODEL  834  A  LIMB  PROFILES  ALONG  ORBIT 


# 
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Figure  2  shows  a  vertical  perspective  projection  of  a  4:30  pm  DMSP  orbit  (this  is  the 
local  time  at  that  point  in  the  orbit  when  the  satellite  crosses  the  equator  going  north).  The 
satellite  altitude  is  850  km.  The  boundaries  of  the  continents  are  positioned  for  a  universal 
time  of  00:00.  The  time  of  year  is  equinox  (either  spring  or  fall)  as  indicated  by  the 
placement  of  the  terminator.  The  SZA  along  the  orbit  is  of  particular  interest  since  it  is  one 
of  the  parameters  defining  the  region  of  useable  data  from  SSULI.  Calculations  of  SZA  have 
been  performed  for  the  above  orbit  and  another  at  summer  solstice.  Figure  3  shows  these 
results  versus  latitude  over  a  full  period  of  the  orbit.  At  equinox,  the  minimum  value  is  about 
67".  The  corresponding  value  at  summer  solstice  is  about  52".  The  range  of  useful  data  is  at 
SZA  values  less  than  approximately  80". 

Simulations  of  SSULI  834  A  limb  observations  (without  noise)  along  the  above 
displayed  orbit  are  given  in  the  lower  panel  of  Figure  4.  The  ionospheric  specification  model 
ICED  (Tascione  et  al.,  1988)  was  used  to  sp«:ify  the  profiles  for  the  multiple  scattering  . 
calculations  (equated  to  ICED’s  electron  density  profile).  The  profiles  corresi>onding  to  the 
834  A  results  appear  in  the  upper  panel.  A  single  model  atmosphere  and  SZA  (70")  have 
been  used  in  order  to  directly  relate  changes  in  834  A  to  O^.  It  should  be  noted  that  these 
results  were  obtained  at  a  later  date  than  those  to  follow  in  subsequent  sections.  The 
radiative  transfer  model  of  Gladstone  (1982)  was  used  to  perform  the  multiple  scattering 
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IHgure  1.  Solar  EUV  spectrum  used  to  calculate  834  A  emission  by  direct  photoionization  and 
by  electron  impact  on  O.  The  spectrum  is  based  on  the  Hinteregger  reference  spectrum  with 
scalings  using  an  F107  value  of  200. 
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Figure  3.  Solar  zenith  angle  along  the  orbital  track  for  4:30  pm  DMSP  orbits  at  summer  solstice 
(upper  panel)  and  at  equinox  (lower  panel). 
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ICEDe  SSN  =  200  EFFQ  =  5  SZA  =  70° 


Hgure  4.  Dayside  ICED  electron  density  profiles  at  selected  latitudes  for  a  4:30  pm  DMSP  orbit 
(upper  panel)  and  corresponding  834  A  limb  profiles  observed  fiom  850  km  (lower  panel). 
The  O*  density  profiles  used  in  the  calculations  were  assumed  to  be  spherically  symmetric  and 
equal  to  the  profiles  in  the  upper  panel.  An  effective  sunspot  number  of  200 ^gnd  an  effective 
Q  value  of  S  were  used  to  q)ecify  the  density  profiles.  A  single  MSIS-86  model  atmo^here  and 
solar  zenith  angle  of  70*  were  used  in  order  to  unambiguously  relate  changes  in  834  A  to  changes 
in  O*. 


calcuh^ons.  Those  to  follow  in  later  sections  used  codes  DSDA  and  LIMB.  Om^aiisons 
have  been  made  between  the  modeb  and  show  similar  results  for  834  A. 

The  shapes  of  the  834  A  profiles  in  Hguie  4  are  similar  in  ^ite  of  die  conreqionding 
profiles  having  values  ranging  fiom  0.5  to  2.5  x  10*  cm  ’  aiul  h^  values  ranging 
from  370  to  500  km.  The  magnitude  of  the  834  A  profile  decreases  as  inoeases  since 
attenuation  by  self  absorption  more  than  compensates  for  increases  in  the  source  function  (see 
Section  5  for  the  dependence  of  the  source  function  on  O^.  To  further  examine  the 
dependence  of  shape  of  the  834  A  limb  {vofile  on  CX,  the  profile  at  -TXf  latitude  (solid 
curve)  has  been  scaled  to  produce  an  834  A  profile  similar  in  magnitude  to  that  at  0*.  The 
upper  panel  in  Figure  5  shows  the  unsealed  CT  profile  at  0"  and  the  scaled  profile  at  -20*. 

The  lower  panel  shows  the  corresponding  834  A  profiles  in  counts  based  on  a  calibration 
factor  of  .2  counts/s/R,  90  s/limb  scan,  and  1  s/measurement  While  the  magnitudes  ate 
similar,  their  peaks  do  differ  by  about  50  km  in  tangent  altitude.  This  difference  arises  for 
profiles  having  the  same  scale  heights.  With  differing  scale  heights,  the  difference  could  be 
either  larger  or  smaller.  The  diamonds  show  a  simulation  of  the  measurements  based  on  the 
solid  curve  using  Poisson  statistics.  Note  that  unlike  the  previous  figure,  the  results  are 
shown  on  a  linear  intensity  scale.  This  simulation  suggests  that  it  will  not  be  possible  to 
distinguish  between  O*  profiles  like  those  in  the  upper  panel  which  ate  in  the  range  expected 
for  daytime  conditions.  These  profiles  show  substantial  differences  with  the  ratio  of 
exceeding  a  factor  of  two  and  the  h^’s  differing  by  more  than  100  km.  The  difference  in 
TEC  up  to  800  km  is  more  than  a  factor  of  two. 


4.  CHAPMAN  FUNCTIONS  USED  TO  CHARACTERIZE  O*  DENSITY  PROFILES 

At  the  request  of  R.  McCoy  of  NRL,  calculations  of  the  834  A  limb  profile  were 
performed  for  O'*  distributions  given  by  Chapman  functions.  These  results  will  be  presented 
in  the  next  section.  McCoy  requested  that  this  be  done  for  five  values  of  N,...  seven  values 
of  h„^,  and  a  single  value  for  H^  (we  selected  this  value  to  be  120  km).  The  N^  values  in 
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Figure  5.  O*  profiles  taken  from  Figure  4  with  a  scaling  of  the  -20*  latitude  profile  (upper  panel) 
in  order  to  produce  similar  magnitudes  in  the  corresponding  834  A  profiles  shown  in  the  tower 
panel.  The  diamonds  represent  simulated  SSUU  data  based  on  the  solid  profile. 
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units  of  lO’  cm'^  are  0.5,  1.0,  S.O,  10,  and  SO.  Units  of  10*  cm'^  will  be  used  in  die  remainder 
of  this  report  when  discussing  N^.  A  value  of  SO  is  somewhat  larger  than  expected  at  solar 
maximum  while  the  values  of  .5  and  1  are  more  impropriate  to  nighttime  condititms.  We  will 
take  note  of  this  when  discussing  the  834  A  results  in  the  next  section.  The  h^  values 
uniformly  span  the  range  from  300  to  600  km  in  SO  km  increments.  The  larger  values  ^iply 
to  (T  distributions  in  the  tropical  arcs  which  are  situated  to  either  side  of  the  geomagnetic 
equator.  Figure  6  shows  Chapman  functions  for  the  seven  h_„  values  and  the  single 
value  of  10  (in  the  above  units). 


5.  EXAMPLES  OF  834  A  VOLUME  EMISSION  RATES 

Examples  of  the  source  function  S(z)  (total  volume  emission  rate)  for  the  long 
wavelength  line  of  the  triplet  (at  834  A)  will  follow.  They  were  obtained  from  code  DSDA 
using  a  single  initial  source  function  S,  (initial  emission  from  the  sources  cited  in  Section  1) 
based  on  the  parameters  given  in  Section  2  and  an  SZA  of  75**.  The  large  SZA  reflects  the 
situation  to  be  encountered  during  a  4:30  pm  orbit  Three  flgures  follow  for  values  of 
0.5,  5.0,  and  50.  Figure  7  shows  S  for  =  .5  and  h^  =  300,  450,  and  500  km  along  with 
the  single  So  cited  above.  Figures  8  and  9  show  similar  results  for  the  other  two  values. 
As  expected,  S  departs  increasingly  from  S,  as  and  h.,,,  become  larger.  Corresponding 
limb  proflles  show  much  less  sensitivity  to  due  to  a  compensating  effect  by  the 
transmission  function  for  self  absorption.  This  was  seen  in  Figures  4  and  5  and  will  be  seen 
in  results  to  follow. 


6.  834  A  LIMB  PROFILES  BASED  ON  CHAPMAN  FUNCTIONS 
Figures  10  •  12  give  834  A  limb  profiles  for  the  combined  three  lines  in  displays 
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=  1 X 10*  cm*^ 


10^  10® 
0+  DENSITY  (cm*^) 


Figure  6.  A  subset  of  the  Oilman  functions  used  to  characterize  O*  for  detennining  the  834 
A  limb  profiles  in  Section  6.  The  values  of  the  Chapman  parameters  for  the  displayed  profiles 
are  =  120  km,  =  10‘  cm  ’,  and  h^  =  300  to  600  km  in  increments  of  50  km. 
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VOLUME  EMISSION  RATE  (cm  *  tec  ') 


VOLUME  EMISSION  RAIE  (cm'*  wc') 


VOLUME  EMISSION  RATE  (cm'*  lec'*) 


Rgure  7.  Source  functions  S  (labeled  as  volume  emission  rates)  for  the  long  wavelength  line  of 
the  834  A  triplet  based  on  O*  profiles  with  =  5x10^  cm*’  and  h^  values  of  300,  450,  and 
600  km.  The  same  initial  source  function  S,  for  an  SZA  of  75"  was  used  calculate  each  of  the 
Sj. 
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Figure  8.  Similar  to  Figure  6  except  for  N„  =  5x10*  cm  *. 
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Figure  9.  Similar  to  Figure  6  except  for  =  5xl0‘  cm"’. 
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Figure  10.  834  A  limb  profiles  versus  look  angle  for  the  three  lines  combined.  The 

corresponding  tangent  altitude  ranges  from  0  to  740  km.  The  profiles  in  the  upper  panel  are  for 
O*  profiles  with  N_.  =  5x10^  cm*^  and  h^  values  fiom  300  to  600  km  in  50  km  increments. 
Similar  profiles  are  displayed  in  the  lower  panel  for  =  10*  cm  *. 
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Figure  11.  Similar  to  Figure  9  except  for  values  of  5x10*  and  10‘  cm  *. 
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convenient  for  showing  the  effect  of  changing  h^.  For  example,  the  upper  panel  of  Figure 
10  shows  profiles  for  the  seven  chosen  values  of  h^  at  the  smallest  selected  value  (0.5). 
The  profiles  are  displayed  versus  look  angle  (0*  refers  to  zenith  viewing)  and  correspond  to 
tangent  altitudes  from  0  to  740  km.  The  remaining  panels  in  these  Hgures  refer  to  the  rest  of 
the  selected  values.  Figures  13  -  16  show  the  same  results  from  the  previous  three 
frgures  expect  versus  in  each  of  seven  panels  for  the  seven  h^  values.  In  a  given  panel, 
the  834  A  profile  is  seen  to  flatten  as  increases.  It  becomes  less  bright  near  and  below 
its  peak  and  brighter  at  high  tangent  altitudes.  The  overall  behavior  is  caused  by  a 
combination  of  changes  in  S  (see  Figures  7-9)  and  transmission  effects  along  the  line-of- 
sight  The  latter  effects  lead  to  the  decrease  in  biighmess  at  the  lower  tangent  altitudes.  At 
high  tangent  altitudes,  these  effects  ate  weak  allowing  for  detection  of  the  enhancement  in  S 
caused  by  increases  in  N^. 


For  remote  sensing  of  0^  measurement  of  changes  in  the  shape  of  the  834  A  is  more 
important  than  measurement  of  absolute  brightness.  Measured  brighmess  is  a  function  of 
calibration,  the  solar  EUV  flux,  neutral  composition,  SZA  along  the  line-of-sight,  and  finally 
O^.  The  uncertainties  associated  with  the  first  three  of  these  items  will  not  allow  for  the 
desired  accuracy  in  determining  (T  from  the  overall  magnitude  of  the  observed  profile. 
Regarding  shape,  there  is  considerable  variadon  over  the  full  range  of  values  considered. 
Nevertheless,  the  variadon  is  weak  over  the  range  expected  for  nominal  daytime  conditions 
(10  to  30). 


7.  CONCLUSIONS 

A  sensitivity  study  has  been  undertaken  to  better  understand  the  expected  properties  of 
834  A  data  to  be  obtained  by  SSULI.  Calculated  834  A  limb  profiles  have  been  presented  for 
O^  i  /files  given  by  the  model  ICED  and  by  Oiapman  functions.  While  considerable 
variation  in  the  shape  of  the  834  A  profile  is  exhibited  over  the  full  range  of  N^  (0.5  to  50 
in  units  of  10*  cm  *)  and  h„,„  (300  to  600  km)  considered  (at  an  SZA  of  75"),  only  a  modest 
change  occurs  over  the  expected  range  of  these  parameters  for  most  daytime  conditions. 
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Figure  13.  Similar  to  Figures  9-11  except  arranged  to  conveniently  compare  profiles  for 
different  values.  Profiles  in  the  upper  and  lower  panel  are  for  values  of  300  and  350 
km,  respectively. 
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Figure  14.  Similar  to  Figure  12  except  for  values  of  4(X)  and  450  km. 


-  20  - 


LOOK  ANGLE  (deg)  LOOK  ANGLE  (deg) 


ft 


INTENSITY  (R) 


100 
102 

104 
106 
108 
110 
112 
114 
116 

118 
120 

100  1000  10000 

INTENSITY  (R) 


Figure  15.  Similar  to  Figure  12  except  for  values  of  500  and  550  km. 
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Figure  16.  Similar  to  Figure  12  excxpt  for  an  value  of  600  km. 
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DEVELOPMENT  OF  AN  INVERSION  ALGORITHM  FOR 
RETRIEVING  MESOSPHERIC  O,  AND  THERMOSPHERIC  Oj 
DENSITY  PROFILES  FROM  SUSIM  SOLAR  OCCULTATION  DATA 
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SUSIM 


CPI  has  been  involved  in  an  ongoing  effort  to  develop  inversion  algorithms  for  the 
retrieval  of  Oj  and  O3  number  density  profiles  fiom  solar  occultation  data  to  be  obtained  from 
the  SUSIM  instrument  on  UARS.  Ttus  problem  is  challenging  because  SUSIM  views  the  full 
solar  disk,  in  contrast  to  most  current  occultation  experiments  in  which  the  instmment  field  of 
view  is  restricted  to  only  a  small  fiaction  of  the  disk.  The  resulting  inversion  problem  requires 
the  solution  of  two  successive  integral  equations,  one  arising  from  the  integration  over  the  finite 
source  and  the  second  due  to  the  integration  along  the  slant  path  through  the  atmo^here,  which 
is  an  inherent  feature  of  the  limb  viewing  geometry.  The  latter  gives  rise  to  the  well  known  Abel 
integral  equation. 

The  initial  method  proposed  for  the  inversion  involved  using  simulations  to  generate  a 
model  dependent  functional  relationship  between  the  extinction  (or  normalized  intensity)  profile 
and  the  slant  column  density  of  the  absorbing  species.  Given  an  extinction  profile  measured  by 
SUSIM,  the  corresponding  column  density  was  obtained  by  a  simple  interpolation  on  this 
functional  relationship,  assuming  an  initial  model  profile  and  iterating  until  some  convergence 
criterion  was  satisfied.  The  number  density  profile  was  then  determined  by  numerically 
evaluating  the  Abel  integral,  using  a  quadratic  least  squares  fit  to  the  column  di..j$ity  in  order  to 
minimize  magnification  of  errors  due  to  statistical  fluctuations  in  the  data.  This  method  was 
capable  of  determining  the  overall  magnitude  of  the  ozone  density  profile  assuming  that  the 
vertical  structure  was  (somehow)  known  beforehand.  This  is  not  usually  the  case  of  course,  and 
in  general  the  algorithm  failed  to  converge  to  the  correct  density  profile  in  realistic  simulations. 


An  improved  inversion  algorithm  has  been  developed  over  the  past  2-3  years  which  shows 
real  promise  for  extracting  useful  Oj  and  0,  vertical  profiles.  This  algorithm  solves  the  finite 
source  integral  equation  by  expanding  both  the  measured  absorption  profile  and  the  unknown 
number  density  profile  in  terms  of  Laguerre  polynomials.  The  integral  equation  is  therefore 
solved  exactly,  resulting  in  a  smooth,  analytic  representation  for  the  density  profile.  This  method 
does  not  require  the  assumption  of  an  initial  model  density  -  the  data  themselves  determine 
unambiguously  the  number  density  profile.  The  essential  assumptions  made  are  that  there  is  only 
one  dominant  absorbing  constituent  at  the  wavelength  of  interest  and  that  the  concentration  of 
this  species  possesses  spherical  symmetry,  i.e.,  it  is  a  function  of  altitude  only. 

A  detailed  discussion  of  the  algorithm  may  be  found  in  the  published  paper  which  has 
been  included  here  in  Appendix  A.  It  is  worthwhile  to  discuss  briefly  several  improvements 
which  have  been  made  in  the  algorithm  since  its  initial  publication.  A  more  general  approach 
has  been  developed  for  the  solution  of  arbitrary  integral  equations  utilizing  the  same  idea  of 
expansions  in  complete  sets  of  orthogonal  polynomials.  This  generalized  algorithm  is  simpler 
both  in  terms  of  the  mathematical  formalism  and  its  numerical  implementation  on  the  computer. 
One  major  practical  improvement  is  the  ease  with  which  solar  limb  darkening  effects  may  be 
incorporated  into  the  new  algorithm.  Other  corrections  to  the  original  solution,  such  as 


wavelength  dependence  of  the  absorption  cross  sections  across  the  instrument  bandpass  and 
variations  of  the  apparent  solar  disk  size  with  altitude  are  also  more  easily  incorporated  into  the 
new  framework. 

By  way  of  example  we  present  in  Figs.  1-4  several  representative  results  from  simulations 
of  SUSIM  occultation  data  for  both  Oj  and  Oj.  These  results  were  generated  using  the  code 
SOLXX!,  which  contains  both  a  forward  module  for  accurate  simulation  of  SUSIM  data  as  well 
as  the  inversion  algorithm  discussed  above.  These  simulations  involve  the  following  steps.  First 
a  model  Oj  or  Oj  density  profile  is  used  to  generate  a  slant  column  density  as  a  function  of 
tangent  altitude.  The  resulting  transmission  profile  is  then  integrated  over  the  solar  disk,  as  well 
as  the  instrument  bandpass,  to  generate  the  integrated  extinction  profile  as  a  function  of  altitude. 
Finally,  a  random  number  generator  is  used  to  simulate  realistic  effects  of  statistical  scatter  in 
the  data.  This  data  profile  is  then  read  into  the  inversion  algorithm  and  the  retrieved  number 
density  profile  is  compared  to  the  original  model  profile  to  determine  the  accuracy  of  the 
inversion  procedure.  It  should  be  noted  that  the  inversion  cannot  be  exact  once  it  has  been 
corrupted  by  random  noise,  which  means  that  it  is  impossible  in  principle  to  extract  exactly  the 
model  density  profile  under  these  conditions.  The  most  one  can  hope  for  is  that  the  model 
profile  and  the  extracted  profile  are  consistent  to  within  the  total  error  bars  specified  on  the 
extracted  profile.  At  this  point  the  error  analysis  has  not  been  completed  and  so  there  are  no 
error  bars  given  on  the  inversion  results. 

In  Fig.  1,  the  top  panel  shows  a  simulated  full  disk  absorption  profile  for  ozone  (open 
squares)  and  the  fit  generated  with  25  terms  in  the  Laguerre  expansion  (solid  line).  The 
wavelength  is  3090  Angstroms,  with  narrow/narrow  slit  configuration.  The  lower  panel  shows 
the  percentage  difference  between  the  data  and  the  fit  as  a  function  of  tangent  altitude.  Note  the 
random  nature  of  the  scatter  about  the  fit,  on  the  order  of  10%  at  high  altitudes.  Figure  2  shows 
the  retrieved  ozone  density  (solid  line)  vs  the  initial  model  density  used  to  generate  the  data 
(dotted  line).  Figures  3  and  4  are  similar  results  for  an  measurement  at  1600  Angstroms. 
Note  that  although  there  is  more  statistical  noise  in  the  Oj  data,  due  to  the  diminished  solar  flux 
at  1600  Angstroms,  the  density  retrieval  is  better  than  for  the  O3  case,  largely  because  the  larger 
scale  height  of  Oj  requires  many  fewer  terms  in  the  Laguerre  expansion  (5  vs  25).  Overall, 
however,  both  results  are  encouraging. 
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HGURE  CAPTIONS 


Hg.  la  Simulation  of  measured  full  disk  absorption  profile  for  ozone  at 
3090  A  (open  squares)  compared  to  the  Laguerre  expansion  fit  with 
N  =  25  (solid  line). 

Fig.  lb  Percentage  difference  between  the  data  and  fit  shown  in  Figure  1. 


Fig.  2  Deduced  ozone  density  obtained  from  retrieval  algorithm  (solid 

line)  compared  to  model  ("true*')  density  used  to  generate  the  data 
(dotted  line). 


Figs.  3a-b  Same  as  Figure  1  except  the  simulated  measurement  is  at  a 
wavelength  of  1600  A  and  the  Laguerre  expansion  has  5  terms. 

Fig.  4  Same  as  Figure  2  but  for  the  Oj  retrieval  problem. 
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AbHtract— An  algorithm  is  presented  for  the  retrieval  of  atmospheric  number  denaty  profiles 
by  inversion  of  full  disk  solar  ooculution  dau.  A  general  derivation  of  the  normalized 
absorption  profile  obtained  from  such  measurements  is  first  given,  assuming  an  arbitrary  solar 
intensity  distribution.  Using  Fourier  transform  techniques,  the  resulting  integral  equation  is 
then  solved  exactly  to  yield  an  expression  for  the  number  density  in  terms  trf  the  measured 
absorption  profile.  This  exact  mathematical  solution  of  the  intepr^  equation  eliminates  the 
need  for  an  iterative  solution  based  on  an  assumed  model  dei^ty  profile  and  makes  the 
inversion  algorithm  presented  here  unique  among  the  currently  used  algorithms  known  to 
the  authors.  The  formal  solution  is  impkmented  numerkally  by  int^ucing  a  suitable 
represenution  of  the  dau  in  temu  of  Laguerre  polynomials,  and  model  cak^tions  are 
presented  for  the  qpedal  case  of  a  uniform  solar  intensity  distribution.  The  results  of  these 
calculations  are  used  to  demonstrate  the  convergence  properties  of  the  algorithm  and  its 
suitability  for  accurate  reproduction  of  several  represenutive  model  density  profiles  with  both 
constant  and  altitude-dqiendent  scale  heights.  Finally,  a  ample  numerical  calculation  is 
performed  to  estimate  the  effect  of  nonunifotm  solar  intensity  distributions  on  the  retrieved 
density. 


1.  INTRODUCTION 

In  this  paper,  we  describe  mathematical  methods  which  can  be  used  to  deduce  the  number  of 
density  profiles  of  atmo^rheric  constituents  from  solar  occultation  measurements.  This  technique 
uses  an  Earth-orbiting  detector  to  measure  the  relative  change  in  solar  irradiance  as  the  Sun  is 
occulted  by  the  limb  of  the  Earth  (either  sunrise  or  sunset).  At  a  given  wavelength,  the  signal  is 
attenuated  due  to  photo-absorption  by  the  species  present  in  the  atmosphere  which  absorb  solar 
radiation  at  that  wavelength.  The  measured  transmission  profile  therefore  contains  information 
about  the  combined  slant  path  column  density  of  all  of  these  species.  In  the  special  case  in  wfaidi 
there  is  one  dominant  absorber  at  the  wavelength  of  interest,  the  data  are  directly  related  to  the 
column  density  of  that  species.  Because  the  occultation  technique  involves  a  relative  measurement, 
with  normalization  given  by  the  measured  solar  irradiance  above  the  atmoq>here  during  eadi  event, 
the  results  are  independent  of  instrument  calibration.  This  is  a  very  advantageous  feature  because 
it  tends  to  negate  the  effects  of  instrument  degradation  over  time.  On  the  other  hand,  the  inversion 
of  full  disk  occultation  data  presents  special  problems  because  the  apparent  size  of  the  source,  as 
seen  by  the  detector,  is  generally  large  compared  to  the  local  atmospheric  scale  height.  A  angle 
measurement  therefore  consists  of  light  received  from  widely  separated  ray  paths,  which  collectively 
sample  a  region  of  the  atmosphere  over  which  the  density  distribution  can  vary  significantly. 

A  good  deal  of  work  performed  in  the  past  was  related  to  stellar  occultation  measurements,'  or 
to  measurements  which  view  only  a  small  fraction  of  the  solar  disk.^  The  calculations  in  this  case 
are  relatively  straightforward  since  a  simple  exponential  relationship  exists  between  the  measured 
transmission  profile  and  the  slant  column  density  at  each  altitude.  In  general,  however,  the  solar 
occultation  problem  requires  an  integration  over  the  finite  source  to  generate  the  transmission 
profile,  and  therefore  the  inversion  process  involves  the  solution  of  a  nonlinear  integral  equation. 
Previous  attempts  at  deducing  species  concentrations  from  full  disk  solar  occultation  data  relied 
on  an  iterative  approach  based  on  the  assumption  of  an  initial  model  density  profile.^  The  integral 


fTo  whom  all  correspondence  should  be  addressed. 
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Here,  I{X,  0,  is  the  unattenuated  solar  flux  (photons/ A/sr/sec)  at  wavelength  I  from  the  source 
point  P  =>  (0,  ^);  9(A)  is  the  absorption  cross  section,  Ntr,(0,  ^)]  is  the  slant  column  density 
corresponding  to  an  integration  along  the  Unenif-sight  from  detector  to  the  point  P,  defined  1^ 
the  tangent  radius  r„  and  Q(X)  is  the  relative  responsivity  fimction  of  the  instrument.  The  geometric 
factor  cos(0  +  C)  gives  the  projection  of  a  light  ray  ori^nating  at  point  P  along  the  line  of  sight 
to  the  detector.  It  should  be  noted  that  f  is  a  function  of  two  variabks:  the  tangent  radius  defined 
by  a  look  angle  V  to  the  center  of  the  source  and  the  observation  wavelength  A. 

The  first  integration  in  Eq.  (1)  is  over  the  bandwidth  of  the  instrument.  The  integration  over 
angles  (0,  is  a  consequence  of  the  finite  source  size  and  represents  the  essential  difference  between 
the  solar  a^  stellar  oocultation  problems.  This  added  complexity  makes  the  inversion  of  the 
extended  source  problem,  i.e.,  the  solution  of  Eq.  (1)  for  N  in  terms  of  F,  much  more  diflBcuh  than 
for  the  point  source  problem. 

Equation  (1)  can  be  simplified  by  making  two  mathematical  approximations  based  on  the 
geometry  of  Fig.  1.  We  first  note  that  K  |  4  I,  where  C».  is  the  solar  half  angle  (~  ISO,  sod 

therefore  the  projection  factor  cos(0  +  C)  %  cos  0.  Second,  we  use  the  fact  that  So>a>r^,  where 
is  the  Sun-to-detector  distance,  a  is  the  solar  radius  and  the  orbit  radius,  to  make  the  approxi¬ 
mation  0^  s  sin~'  (1  —  (a/So'fY’^  st  x/2.  Finally,  we  make  the  physical  assumption  of  azimuthal 
synunetry  in  the  solar  flux  distribution,  so  that  /(A,  0,<^)  = /(A,0).  Equation  (1)  then  becomes 

f(A.  =  dAe(A)  f  d0  sin  0  cos  0/(A,  0)  j  ^  (2) 

F  j/u  Jo  Jo 


and 


K  =  2nj  dAe(A)  d0  sin  0  cos  0/(A,  0), 

where  we  have  defined  the  quantity  =r,a/So.  with  r,  =  (rj— ri)*'*  (see  Fig.  1). 

This  result  may  be  rewritten  in  a  more  transparent  form  by  introducing  a  change  of  variables 
from  (0,  (ft)  to  (ji—  cos  0,  y  =  0  sin  0  cos  0)  and  interchanging  the  order  of  integrations  to  obtain 


F{z) 


|^dAe(A) 


(3) 
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the  meMurement  region  extends  to  much  higher  altitudes,  it »  straightforward  to  include  the  next 
tom  in  the  expansion,  of  order  z/r^,  in  Eq.  (9).  Finally,  it  is  convenient  to  make  one  more  diange 
of  variables  to  the  dimensionless  quantity  m  » [(r '  —  z  which  leads  to  the  following  expression 
for  the  number  density: 


"(*)- 


-1 
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(10) 


For  numerical  computation,  Eq.  (10)  has  an  advantage  over  Eq.  (9)  in  that  there  is  no  singularity 
at  the  lower  limit. 

In  concluding  this  section,  we  wish  to  point  out  that  Eqs.  (3)-(S)  can  also  be  obtained  if  the 
radiation  source  is  treated  as  a  two-dimensional  disk  ratto  than  a  qtbere,  as  was  done  here. 
The  equivalenoe  of  the  two  geometries  can  be  established  provided  that  C.  <  1>  t-e.,  the  ratio  of 
the  nu^tion  source  size  to  the  observing  distance  is  mudi  lea  than  unity.  For  sat^te  observations 
of  the  Sun,  C«  4  x  10~^,  and  therefore  the  solar  disk  geometry  will  be  an  excdlent  approximation 
to  the  exact  treatment  given  here.  Roble  and  Norton^  used  the  solar  disk  getHnetiy  and  obtained 
an  expression  similar  to  Eq.  (3).  However,  the  function  /f(i.y)  defined  in  Eq.  (4)  appears  to  differ 
from  their  result  except  in  the  limit  of  a  uniform  solar  intensity  distribution,  /(A,ii)-»/(A).  This 
difiereiKe  results  from  different  interpretations  concerning  the  luture  of  the  intensity  distribution 
acroa  the  Sun.  We  believe  that  the  derivation  presented  here  is  the  most  general  possible  and  is 
consistent  with  actual  solar  observations. 


^  3.  FOURIER  TRANSFORM  SOLUTION 

For  the  purpose  of  the  present  work,  we  shall  make  the  assumption  that  the  variation  of  e(A) 
over  the  instrument  bandwidth  in  Eq.  (3)  can  be  ignored,  or  suitably  averaged,  in  performing  the 
inversion.  The  validity  of  this  assumption  depends  on  the  interplay  between  two  factors:  the  exact 
bandwidth  of  the  instrument  and  tiw  detailed  structure  of  the  absorption  cross  section  over  this 
bandwidth.  In  practice  this  must  be  evaluated  for  the  particular  case  of  interest.  Assuming  a 
9  sufiSdently  narrow  bandwidth  and  a  smoothly  varying  <r(A),  however,  the  approximation  is 
reasonal^.  While  this  assumption  simplifies  the  solution  of  Eq.  (3)  significantly,  it  is  not  essential 
to  the  inversion  algorithm  presented  in  the  following  sections.  A  straightforward  generalization  of 
the  following  solution  which  allows  for  the  full  functional  dependence  of  e(A)  will  be  given 
elsewhere.*  For  the  present,  we  will  assume  that  <r(A)  =  a,,  a  constant  evaluated  at  the  center  of 
the  bandpass.  The  function  1(1,  n)  can  then  also  be  njptaced  in  Eqs.  (4)  aixl  (5)  by  its  value  at  the 
central  wavelength,  or  one  may  use  the  fully  folded  solar  spectrum  defined  by 

/4i(A<)  =  j^<WG(-()/(^M).  (11) 

For  simplicity  of  notation,  we  will  subsequently  denote  this  effective  solar  intensity  distribution  by 
Note  that  with  this  assumption  made,  the  explicit  wavelength  integration  no  longer  appears 
in  Eqs.  (3)  and  (S)  and  the  kenwl  is  now  a  function  of  the  variable  y  alone,  i.e.,  mi.,y)-*H{y). 
i  In  order  to  facilitate  the  mathematical  development,  we  next  introduce  the  absorption  profile 

defined  by 

C(z)=l-F(z).  (12) 

This  function  has  the  useful  property  that  C(oo)  =  0  since  F(oo)  =  1.  Equation  (3)  is  then  rewritten 
in  the  following  form: 

^  G(z)  =  j*^dyH(y}x(z+y).  (13) 


where 


» 


X(u)  =  I  —  e"**'*'** 


(14) 
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From  Eq.  (14).  the  column  density  is  given  by 

iV(z)--ilnll-y(z)l 

and,  therefore. 

<W(z)  ^  1  1  dy(z) 

d*  ffo  I  —  X(*)  dz 


Inserting  this  result  mto  Eq.  (10),  the  expression  for  the  number  density  becomes 


«(z)  = 


(24) 

(25) 

(26) 


where  the  prime  indicates  differentiation  with  respect  to  the  complete  argument  of  y.  Since  y  has 
been  expressed  in  terms  of  the  Fourier  transform  of  the  measured  extinction  profile,  Eq.  (26) 
constitutes  the  desired  solution  for  the  number  density  as  a  function  of  tangent  altitude. 


4.  EXPANSION  IN  LAGUERRE  POLYNOMIALS 
It  is  obvious  from  the  results  of  the  last  section,  and  in  particular  from  Eq.  (26).  that  the  retrieval 
of  accurate  niunber  density  profiles  from  solar  occultation  measurements  requires  an  accurate 
evaluation  of  the  integral  in  (23).  Although  this  integral  looks  relatively  simple,  it  is  in  general 
highly  singular.  To  see  this,  we  consider  the  following  polynomial  parameterization  of  the  solar 
limb  darkening  function: 


/(A,M)//(A,l)-iay(A)/i^,  (27) 

y-» 

where  the  coefficients  aj{X)  are  empirically  determined.  Consistent  with  the  discussion  in  Sec.  3, 
these  will  be  replaced  by  their  average  value  over  instrument  bandwidth,  and  denoted  simply  by 
Oj.  The  kernel  then  takes  the  form 

').  (28) 

where  B{x,y)  is  the  beta  function.  The  corresponding  futKtion  Y{p)  becomes 


(29) 


Here,  J,  is  the  Bessel  futKtion  of  order  v,  r(x)  is  the  gamma  fuiKtion  and  the  normalization  factor 
is  *  n'LjajHj  +  2). 

The  function  Y(j>)  has,  in  general,  an  infinite  number  of  real  zeros  pi,  which  means  that  the 
integrand  in  Eq.  (23)  has  an  infinite  number  of  simple  poles  on  the  real  axis.  In  addition  to  the 
there  will  also  be  singularities  in  the  function  Cl(p)  and  the  integral  in  Eq.  (23)  is  determitKd 
by  the  sum  of  residues  at  these  poles.  Schematically,  we  write  this  sum  as  follows: 


Xiz)  =  2ni^  + 1  =  y‘'>(z)  +  y«>(z). 


(30) 


where  Rj  and  R^  are  the  residues  evaluated  at  the  poles  of  S{p)  and  \/Y{p),  respectively.  Since 
z  ^  0  by  definition,  the  contour  is  closed  in  the  upper  half  plane  and  only  poles  on  the  real  axis 
or  with  positive  imaginary  parts  contribute  to  the  integral.  The  contribution  from  the  poles  of 
\IY{p)  is  of  the  form 


y<'‘(r)  =  C5:e^ 

t 


,0iPk) 

Y\p,y 


(31) 
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where  ^.(p)  is  the  Fourier  transform  of  V,{z).  This  sum  may  be  evaluated  analytically  using 

(i-ir 


r 


d«e-L,(u)  ■ 


with  the  result 


i  iip-\r 


(39) 


(40) 


The  function  C(p)  therefore  has  a  single  pole  of  order  n  +  1  at  p  =  il2.  The  reader  should  keep 
in  mind  that  the  conjugate  variable  p  is  also  now  a  dimensionl^  quantity,  consistent  with  the 
substitution  p  -*  fip.  Combining  Eqs.  (23),  (38)  and  (40),  we  can  write  y(z)  in  the  form 


where  the  function  r,(/)  is  defined  by 

r.(/) 


J_  f" 

2«J-. 


dpe^ 


1  (»p-lr 

Y(p)W+Tr~' 


In  order  to  evaluate  this  integral,  we  write 


cu) 


(41) 


(42) 


(43) 


«•  1/2 


and,  therefore. 


2n  nlVdtoj  J.«  ^^y(p)  /p  + 

nlVdo^  Yi 


-ir 

a> 


■1/2 


(«+ir 


(44) 


The  integral  in  Eq.  (44)  has  been  evaluated  by  taking  the  residue  at  the  single  pole  p  =  ko,  in 
accordance  with  the  preceding  discussion. 

The  function  r,(t)  can  be  evaluated  directly  from  Eq.  (44),  wluch  results  in  a  closed  form  solution 
in  terms  of  associate  Laguerre  polynomials,  Z.^’(r).  However,  a  more  consistent  approach  is  to 
expand  r,(/)  in  terms  of  Laguerre  functions,  as  we  did  for  G{t).  We  therefore  write 


r.(0=  I  yiv.it). 


(45) 


with 


y;  =  j’drr.(/)!P*(/).  (46) 

Inserting  Eq.  (44)  into  Eq.  (46)  and  using  Eq.  (39)  to  evaluate  the  integral,  pves  the  following  result: 

(47) 


.-1/2 


It  is  obvious  that  yj  =  0  for  A:  ^  n,  since  the  factor  (<u  -  1/2)  then  gives  zero  in  the  limit  a  -*  1/2. 
The  expansion  in  Eq.  (45)  is  therefore  terminated  at  k  =  n.  Upon  further  evaluation,  it  can  be 
shown  that  the  coefficients  yj  depend  only  on  the  difference  («  —  k)  and  so  can  be  written  as  a 
function  of  one  index  only.  The  reader  is  refened  to  Appendix  1  for  a  detailed  derivation  of  these 
coefficients.  The  results  can  be  summarized  as  follows: 

r.(/)=i:  y/^,.//), 

/-o 


(48) 
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expansion  coefficients  and  the  functions  r.(/)  and  r;(i)  (the  prime  indicates  differentiation)  are 
then  calculated  to  order  M  using  Eqs.  (37X  (48)  and  (S3),  respectively.  Once  these  quantities 
are  known,  the  functions  y(/)  and  genoated  from  the  expansions  (41)  and  (SI)  and.  finally, 

the  Abel  integral  (Eq.  (26)]  is  evaluated  numerically  to  obtain  the  number  density  profile.  The  data 
enter  only  through  the  coefficients  g.,  and  this  is  therefore  the  only  part  of  the  calculation  which 
is  specific  to  a  particular  measurement.  The  /'.(r)  and  their  derivatives  are  characteristic  functions 
of  the  problem;  they  are  completely  determined  by  the  kernel  of  the  integral  equation  (13)  and  thus, 
ultimately,  by  the  form  of  the  solar  intensity  distribution  through  Eq.  (4).  It  is  only  necessary  to 
calculate  them  once  for  a  given  parameterization  of  /^(/i),  i.e.,  for  a  given  wavelength.  This 
calculation  can  be  done  as  precisely  as  possible  and  the  resulu  written  to  a  file  for  storage.  In  the 
actual  daU  analysis,  one  need  only  calculate  the  g.  to  the  desired  order,  read  in  r,(/)  and  r‘,(t) 
to  the  same  order  and  then  do  the  Abel  integration.  This  is  a  very  desirable  feature  of  the  algorithm 
in  terms  of  accuracy  and  computational  efficiency. 

The  calculation  of  the  expansion  coeffidenU  g.  is  a  critical  part  of  the  inversion  algorithm,  and 
it  is  therefore  worthwhile  to  briefly  summarize  the  method  used  to  evaluate  the  integral  in  Eq.  (37). 
Because  the  data  only  determine  (7 (r)  up  to  a  maximum  altitude,  some  physical  assumption 
must  be  made  regarding  the  behavior  of  this  function  as  r  -» oo.  We  assume  an  exponential  falloff, 
so  that  C(t)  =  C(f^)exp(— a(t  — 1,„)]  for  t  '^t^.  The  parameter  a  may  be  determined  by 
actually  fitting  the  data  to  this  form  at  or  it  may,  for  example,  be  assumed  to  be  given  by 
CL  where  H  is  the  (assumed)  scale  height  of  the  absorbing  species.  In  practice,  the  resulte 

show  little  sensitivity  to  the  choice  of  a.  The  integral  in  Eq.  (37)  is  evaluated  by  breaking  the 
integration  region  into  intervals  determined  by  successive  zeros  of  V'.(().  and  integrating  over  each 
interval  using  an  adaptive  Romberg  integration  scheme.  In  performing  the  integration  it  is 
necessary  to  interpolate  the  measured  profile  G{t)  from  the  data  grid  onto  the  integration  grid. 
Because  realistic  dau  always  contein  some  statistical  scatter,  or  noise,  the  profile  G(t)  will  not  in 
general  be  smooth  in  the  sense  of  an  aiuUytic  function.  It  is  important  to  realize  that,  in  the  limit 
as  the  number  of  terms  in  Eq.  (34)  goes  to  infinity,  this  noise  in  the  dau  will  be  reproduced  exactly. 
By  truncating  the  expansion  to  a  smaller  number  of  terms,  we  are  effectively  low  pass  filtering  the 
data,  which  allows  us  to  control  the  amount  of  noise  which  is  passed  on  to  the  density  profile. 
Treatment  of  noisy  data  and  a  deUiled  discussion  of  the  associated  error  analysis  will  not  be 
itKiuded  here  but  will  be  presented  in  a  forthcoming  paper.*  In  the  present  calculations  the  function 
G(/)  is  smooth  by  construction  (see  discussion  in  S.2)  and  it  can  be  accurately  evaluated  at 
each  integration  point  by  doing  a  quadratic  interpolation  on  the  log  of  the  daU  profile.  As  the 
results  presented  below  demonstrate,  the  above  procedure  yields  a  very  accurate  determination  of 
the  g,. 

For  the  purpose  of  illustration,  we  will  consider  the  case  of  a  uniform  solar  intensity  distribution 
in  the  results  adiicb  follow.  While  this  is  the  simplest  possible  case,  it  contains  all  the  eMenrul 
features  of  the  solution  and  is  sufficient  to  demonstrate  the  convergence  properties  and  the  accuracy 
of  the  calculation.  For  this  case  the  function  Y(p)  =  2y,(p)/p  and  therefore  the  zeros  needed  to 
calculate  the  coefficients  jj  and  Cj  [see  Eqs.  (49)  and  (S4c)],  arc  simply  the  zeros  of  the  first  order 
Bessel  function,  vdiicb  are  known  analytically.  In  general,  the  roots  of  the  equation  Y(p)  =  0  must 
be  determined  numerically  siiue  no  analytic  solution  exists  for  arbitrary  T(p)  as  given  by  Eq.  (29). 
Note  that  in  the  present  case  Y’{p„)  is  proportional  to  /o(P*)  =  (-  l)*l-fo(P*)l.  so  that  the  series 
in  Eqs.  (49)  and  (S4c)  altentate.  Convergence  is  therefore  slow  and  it  is  essential  that  the  summation 
be  done  exactly  to  all  orders  in  order  to  obtain  accurate  results  for  the  coefficients  jj  and  Cj.  In 
Appendix  2  we  briefly  discuss  the  method  used  to  sum  the  series  in  Eq.  (49);  the  same  arguments 
apply  to  the  evaluation  of  Eq.  (S4c). 

In  Figs.  2  and  3  the  functions  F.fr)  and  r'(r)  are  plotted  for  several  representative  values  of 
n.  The  maximum  value  t  —  5  corresponds  to  an  altitude  range  —  z^  =  5p  tsSO  km.  These 
functions  retain  the  general  oscillatory  behavior  characteristic  of  the  y',(r);  the  frequeiKy  of 
oscillation  is  higher  for  larger  n  and,  for  a  given  order,  decreases  as  .'  increases.  However,  unlike 
the  which  have  an  dependence,  the  r,(r)  arc  seen  to  increase  with  increasing  ordcr.“- 
In  addition,  for  a  given  value  of  n,  r.(/)  is  in  general  substantially  larger  than  9',(r),  especially 
for  large  n:  Because  of- the  different  behavior  of  the  two  functions  in  the  limit  n  -» oo,  the  series 
in  Eq.  (41)  will  not  converge  as  fast' as  the  original  expansion  for  the  data  given  by  Eq.  (34),  and 
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Fig.  3.  DerivatiYc  of  the  functioiis  r,(i)  for  the  case  of  a  unifonn  solar  intensity  distribution. 


representative  density  profiles.  A  systematic  application  of  the  algorithm  to  the  analysis 
or  realistically  simulated  data  for  specific  atmospheric  constitwnts  is  deferred  to  a  forthcoming 
paper. 

We  first  consider  the  case  of  a  purely  exponential  model  density  profile  of  the  form 
n^(z)  =  iioexp(-r//0.  where  the  scale  height  N  is  independent  of  altitude.  This  is  the  simplest 
realistic  form  of  the  model  density  and  provides  the  least  stringent  test  of  the  algorithm,  since  there 
is  no  vertical  structure  to  reproduce  in  the  profile.  The  most  important  physical  parameter  in  this 
problem  is  the  ratio  S  =  which  is  a  measure  of  how  much  the  density  varies  at  the  tangent 
point  due  to  the  finite  extent  of  the  Sun.  One  would  expect  that  as  3  -» co.  the  normalized 
absorption  profile  yields  the  spedcs  cohmm  density  directly,  since  C(f)»  ;t(/)  from  Eq.  (13).  In 
this  limit  the  solar  occultation  problem  becomes  equivalent  to  the  point  source  problem. 

In  all  the  examples  to  follow,  the  absorption  cross  section  is  chosen  to  be  Og  =  I  in  order  to  make 
the  calculations  as  general  as  possible.  The  effect  of  this  normalization  is  that  the  column  density 
is  now  expressed  in  units  of  Ug as  is  obvious  from  setting  <ro=  1  in  Eq.  (3).  The  corresponding 
normalized  number  density  is  expressed  in  units  of  (see  Eq.  (26)].  These  calculations 

may  therefore  be  extrapolated  to  the  case  of  absorption  by  a  particular  species  simply  by- 
multiplying  the  resulting  number  density  by  this  numerical  factor,  with  a  cross  section  appropriate 
for  that  species  at  the  wavelength  of  interest  As  a  first  example,  we  choose  a  scale  height  H  =  5  km, 
which  implies  3  as  1/5.  The  normalized  model  density  profile  and  the  resulting  absorption  profile 
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dC(2)/G(x)  (5S) 


Fig.  S.  (a)  Pncentags  error,  (<7ti(r)  -  (r(z)]/G(r),  where  G(z)  is  the  dau  profik  plotted  in  Fig.  4(b)  and 
Gk  (z)  is  obtained  the  acpansion  (34)  wild  Mr -40, 80  a^  120  terms,  (b)  Corresponding  percentage 
errors  in  the  (notmalind)  deduced  number  density  as  obtained  from  Eq.  (2d). 

In  order  to  study  the  convergence  of  the  expansions  in  more  detail,  the  quantities  SG(i)  and  SxU) 
defined  by 

5(?(f )  =  <?(<)- I  (55) 

mm9 

and 

SxO)  =  X-oe(0  -  Z  (56) 

■  aO 

are  shown  in  Table  1  as  a  function  of  M  at  several  altitudes.  Here,  Xmo^iO  is  given  by  Eq.  (14) 
with  N{t)  s  N^(tY  At  a  fixed  altitude,  the  addition  of  successive  terms  in  the  expansion  results 
in  errors  which  oscillate  about  zero  with  decreasing  magnitude.  Note  that  SG{t)  is  still  converging 
to  zero  at  Af  a  120,  whereas  in  general  Sx(t)  reaches  a  minimum  value  at  some  lower  order  in  the 
expansion  and  does  not  change  as  more  terms  are  added.  This  saturation  is  a  consequence  of  the 
precision  of  the  jj  coefficients,  which  is  determined  by  the  choice  of  the  cutoff  parameter  Nj  used- 
to  evaluate  the  sum  in  Eq.  (49)  (see  Appendix  2).  As  Nj  is  increased,  the  asymptotic  forms  used 
to  approxinute  the  sum  become  more  precise.  However,  there  is  always  some  finite  error  in  the 
yj,  rmd  this  translates  into  an  dfective  lower  bound  on  the  error  Sx(.0- 
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dn(z)/n(z)  (35) 


Fig.  6.  (a)  As  in  Fig.  S(a)  except  that  the  data  profile  has  been  generated  using  a  model  density  with  a 
12  lun  scale  height  and  the  fit  is  obtained  from  Bq.  (34)  with  A/  «>  10, 20  and  40  terms,  (b)  Corresponding 
percentage  errors  in  the  dedu^  number  density  profile. 


and  no  was  chosen  to  yield  a  peak  to  minimum  ratio  of  10  at  z  =  Zg.  The  data  profile  generated 
by  this  model  density  has  been  plotted  in  Fig.  8(a).  There  is  a  good  deal  more  structure  evident 
in  this  absorption  profile  than  in  the  previous  example,  although  it  is  still  very  subtle  in  comparison 
to  the  structure  in  the  density  itself  [see  Fig.  9(a)].  The  dotted  curve  in  Fig.  8(a)  represents  the  fit 
to  G(z)  obtained  with  M  =  100,  and  Fig.  8(b)  shows  the  corresponding  percentage  error,  which 
0  is  within  1%  over  most  of  the  altitude  range. 

In  Fig.  9(a),  the  deduced  number  density  profile  is  compared  to  the  model  profile  for  the 
cases  A/  =  100,  150  and  200,  with  corresponding  percentage  errors  plotted  in  Fig.  9(b).  With  100 
terms  in  the  expansion,  the  deduced  density  underpredicts  the  sharp  minimum  at  75  km.  The  fit 
in  this  region  is  much  better  when  the  expansion  order  is  increased  but,  as  Fig.  9(b)  indicates,  even 
with  M  =  200,  the  maximum  error  is  still  on  the  order  of  10%  at  this  point  in  the  profile. 
Nevertheless,  these  results  demonstrate  that  the  inversion  algorithm  is  capable  of  retrieving  more 
complicated  model  density  profiles,  although  in  general  significantly  more  terms  are  required  in  the 
expansion  and  the  accuracy  is  somewhat  diminished.  The  important  point  is  that  the  expansion 
is  still  converging  so  that,  in  principle,  good  accuracy  can  be  achieved  even  in  the  case  of  profiles^ 
with  significant  vertical  structure.  Several  methods  are  currently  being  investigated  which  might  be 
used  to  improve  convergence  of  the  expansions  in  these  cases,  so  that  fewer  terms  will  be 
required. 
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MIDDLE  ATMOSPHERE  CLIMATOLOGY 


Over  the  past  three  years  Mr.  William  Sawchuck  of  CPI  has  provided  theoretical  support 
to  the  Middle  Atmosphere  Program  (MAP)  at  the  Naval  Research  Laboratory.  While  most  of  this 
effort  has  been  directed  toward  the  area  of  middle  atmospheric  research.  He  has  also  worked  on 
several  modeling  projects  to  study  the  atmospheres  of  Titan,  lo,  Triton,  and  Uranus.  Primary 
contributions  are  divided  into  three  general  areas  as  outlined  below. 


Software  Development  (Data  Analysis,  Comparison,  and  Reduction) 

The  MAP  theoretical  effort  at  NRL  is  currently  directed  toward  the  characterization  of 
the  spatial  and  temporal  distributions  and  variability  of  chemical  species,  such  as  ozone  and  water 
vapor,  in  the  terrestrial  stratosphere  and  mesosphere.  These  theoretical  efforts  have  relied  on  the 
application  of  several  state-of-the-art  computer  models  which  are  used  to  simulate  the  latitudinal 
and  seasonal  distributions  of  these  important  trace  gases  in  the  atmosphere.  A  large  component 
of  Mr.  Sawchuck’s  work  at  NRL  has  been  focused  toward  the  development  of  software  which 
is  used  in  conjunction  with  these  theoretical  models  to  analyze  and  interpret  their  output  and  in 
many  cases  compare  them  to  the  results  of  other  models  and  observational  data  obtained  from 
ground-based,  satellite,  and  spacecraft  measurements.  Some  of  this  software  is  also  used  to 
reduce  large  observational  data  sets  for  use  as  model  input  data.  In  addition,  many  specialized 
model  analysis  routines  as  well  as  several  generalized  graphical  interface  routines  have  been 
written  for  use  with  the  MAP  atmospheric  models .  A  typical  example  which  was  produced  with 
one  of  these  generalized  routines  can  be  seen  in  figure  1. 


Software  Development/Maintenance  (Theoretical  Models) 

Over  the  course  of  this  contract  the  MAP  group  has  maintained  several  large  atmospheric 
models.  These  include  a  1-Dimensional  time  dependent  photochemical  model  which  incorporates 
vertical  transport  by  eddy  and  molecular  diffusion,  a  general  2-Dimensional  dynamics  model 
which  calculates  meridional  wind  fields  from  atmospheric  temperature  data,  and  a  1-Dimensional 
neutral  and  ionospheric  planetary  model.  Over  the  last  few  years,  Mr.  Sawchuck  has  been 
involved  with  the  modification  and  continuing  development  of  all  of  these  models. 

One  of  the  primary  tasks  in  this  area  is  the  development  of  a  generalized  2-Dimensiona] 
Chemical/Dynamical  model.  The  approach  to  this  problem  is  to  merge  the  above  mentioned  1-D 
photochemical  model  with  the  2-D  dynamical  model.  This  includes  coding  of  many  of  the 
routines  which  are  related  to  the  incorporation  of  a  self-consistent  treatment  of  chemistry  into 
the  dynamical  model.  As  an  example,  some  of  these  routines,  based  on  the  1-D  model,  have  been 
designed  to  read  an  input  list  of  chemical  reactions  and  reaction  rate  constants  and  then 
calculate  the  effective  production  and  loss  rates  for  each  chemical  species.  Since  the 
2-D  chemical/dynamical  model  is  computationally  intensive,  it  was  necessary  to  invert  the  code 
to  run  on  the  NRL  Cray  X-MP/24. 
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Project  Participation  (Theoretical  Modeling) 

Mr.  Sawchuck  has  participated  in  several  modeling  projects  during  the  past  three  years. 
These  have  included  investigations  on  the  seasonal  variation  of  water  vapor  and  ozone  in  the 
upper  mesosphere  and  the  diurnal  variation  of  mesospheric  ozone  (figure  2).  Also,  he  has  been 
extensively  involved  in  the  generation  and  development  of  a  trace  constituent  climatology. 

Responsibilities  included  all  of  the  model  calculations  made  for  the  seasonal  variation 
project  as  well  as  that  of  the  diurnal  ozone  project.  These  required  the  development  of  the  model 
comparison  and  analysis  software  for  both  of  tl^se  projects.  The  paper  that  resulted  from  the 
H20/03  seasonal  variation  work  won  a  1990  Annual  Research  Publication  Award  at  NRL. 

The  development  and  generation  of  a  reference  climatology  of  trace  constituents  in  the 
atmosphere  is  based  upon  the  best  observational  data  available,  and  upon  model  extrapolations 
of  these  data  where  measurements  are  not  currently  available.  We  are  now  entering  the  third 
year  of  this  project  and  the  database  currently  includes  most  of  the  minor  constituents  in  the 
atmosphere  which  affect  UV,  visible,  and  IR  backgrounds,  i.e.,0,  02,  N2,  NO,  N02,  N20, 
HN03,  CO,  C02,  H20,  and  03  (see  figures  3  &  4).  This  work  has  recently  been  presented  at 
a  conference  of  the  International  Union  of  Geodesy  and  Geophysics  (lUGG)  in  Vienna,  Austria 
in  August  1991.  It  is  expected  that  this  work  will  be  published  when  completed. 
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HGURE  CAPTIONS 


Figure  1.  Results  firom  the  Triton  Atmospheric  Model  which  show  methane  tangential  column 
densities  as  a  function  of  altitude  and  pressure.  T(^  panel  is  for  summer  hemisphere  at  45 
degrees  latitude  and  bottom  panel  is  for  winter  hemi^here  at  26  degrees  latitude.  Solid  lines  ate 
Voyager  2  UVS  measurements.  Other  lines  ate  model  calculations  for  various  values  of  the  eddy 
diffusion  coefficient  K.  (Strobel  and  Summers,  1990) 


Figure  2.  Relative  diurnal  variation  of  the  0.4  to  0.8  MHz  brightness  temperature  difference 
normalized  by  the  average  value  over  1200  +/>  2hr  time  period.  Thick  solid  line  is  ground  based 
microwave  data,  other  lines  represent  1-D  photochemical  model  calculations  based  on  various 
suggested  modifications  to  standard  mesospheric  odd  hydrogen  and  odd  oxygen  photochemistry. 
(Zommerfelds  et  al.,  1989) 


Figure  3.  Flow  chart  illustrating  the  general  process  by  which  the  Trace  Constituent  Qimatology 
Database  was  constructed. 


Figure  4.  Complete  zonally  averaged  climatology  for  ozone.  (Summers  and  Sawchuck,  1990) 
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Fig.  8.  (a)  Solid  curve  icpnsenta  the  data  profile  geoeraled  miiig  the  model  density  pkMied  in  Fig.  9(a). 
Dotted  curve  leptesentt  tile  fit  to  the  data  obtained  from  Eq.  (34)  with  100  terms,  (b)  Percentage  error 
in  the  fit  to  the  data  using  100  terms  in  Eq.  (34). 


where  /^(l}  is  the  intensity  at  the  center  of  the  disk.  This  representation  corresponds  to  taking  m  =  ) 
in  Eq.  (27).  The  corresponding  value  of  obtained  from  Eq.  (29)  after  some  manipulation,  is 


(60) 


where  I,(x)  is  the  first  order  modified  Bessel  function.  Solar  u.v.  limb  darkening  measurements* 
give  a  v^ue  B  =  0.797  at  2600  A,  whidi  implies  an  intensity  at  the  limb  (m  =>  0)  which  is  only  20% 
of  the  value  at  the  center.  In  order  to  determine  the  eflfect  of  this  variation  on  the  deduced  number 
density,  a  data  |m>file  was  generated  using  a  single  term  in  Eq.  (S8),  with  bt  =  1  and  a,  =  2.S.  This 
profile  was  thm  inverted  to  obtain  number  density  profila  both  with  and  without  the  limb 
darkening  effect  included.  The  resulting  density  profiles  are  plotted  in  Fig.  10,  where  the  solid  line 
corresponds  to  taking  B  =0  (no  limb  darkening  included)  and  the  dotted  curve  includes  the  limb 
darkening  effect  by  using  the  value  of  B  quoted  above.  It  is  obvious  that  the  effect  on  the  deduced 
density  is  substantial;  the  two  curves  differ  by  approx.  50%.  The  analyns,  which  ignores  limb 
darkening  effects,  results  in  a  larger  deduced  density  because  it  effectively  sees  more  photons  and,_ 
in  order  to  account  for  the  measured  amount  of  absorption,  most  deduce  a  larger  density.  This 
simple  result  demonstrates  that  it  is  essential  to  account  accurately  for  realistic  solar  intensity 
distributions  in  the  analysis  of  solar  occultation  data.  A  more  systematic  study  of  these  effects, 
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Fig.  10.  Deduced  number  density  profiles  obtained  from  the  exponenttal  model  data  profile  of  Eq.  (38) 
with  y  »  I.  The  dotted  curve  resulu  from  using  a  meamred  a.v.  hmb  darkcamg  coHlirient  B  in  Eq.  (60^ 
whereas  the  solid  curve  is  the  result  obtained  by  ignoring  fimb  dailcening  effects. 

As  an  example  of  the  numerical  calculations,  an  algorithm  has  been  successfully  devdoped  to 
extract  density  profiles  in  the  case  of  a  uniform  solar  intensity  distribution.  Numerical  results 
of  model  calculations  presented  here  demonstrate  that  the  algorithm  is  capable  of  accuratdy 
reproducing  physically  motivated  model  density  profiles  used  to  generate  model  data.  In  general, 
the  expansions  which  represent  the  extracted  column  density  are  found  to  converge  more  rapidly 
in  the  case  where  the  density  profile  does  not  contain  significant  vertical  structure.  However, 
accurate  results  can  be  obtained  even  in  cases  where  there  is  large  structure  if  enough  terms  are 
allowed  in  the  expansion.  There  ate  currently  several  options  being  considered  which  might  be 
implemented  in  an  effort  to  improve  the  convergence  a^  accuracy  of  the  expansions. 

The  computer  code  which  has  been  developed  to  do  the  calculations  for  a  uniform  solar  intensity 
distribution  presented  here  can  be  generalized  to  do  the  analysis  for  arbitrary  intensity  distributions 
outlined  in  Secs.  3  and  4  of  the  paper.  Preliminary  results  from  the  present  analysis  indicate  that 
the  incluaon  of  realistic  limb  darkening  effects  in  the  calculation  has  a  significant  impact  on  the 
deduced  number  density  profile  and  is  therefore  essential  to  the  retrieval  of  accurate  den^ty  profiles 
from  solar  occultation  data.  This  question  obviously  deserves  careful  numerical  study  and  efforts 
are  currently  underway  to  make  the  necessary  modifications  to  the  code. 

In  conclusion,  we  have  developed  a  rigorous  mathematical  method  and  accurate  numerical 
procedures  to  deduce  number  density  profiles  of  atmospheric  spedes  from  solar  occultation 
measurements.  This  aiqyroach  can  be  extended  to  include  the  analysis  of  multichannel  measure¬ 
ments  capable  of  extracting  information  about  several  species  simultaneously.  It  is  also  hoped  that 
it  may  be  used  as  a  starting  point  to  attack  the  problem  of  nonspherical  density  distributions  using 
a  sin^r  measurement. 
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by  symmeiry,  since  p_.  »■  —  p,  and  1"( — p)  =  —  Y\p).  This  follows  directly  from  the  fact  that  H(y) 
defined  by  Eq.  (6)  is  an  even  function  of  y.  Using  these  results,  Eq.  (A9}  becomes 


1  i  ‘  r  *  I  '1 
Y(p)  ,h,Y\p.)lp-p,  P.i 


(All) 


and.  therefore. 


\lFb[iornF7.  V.II-. 

The  first  and  third  terms  of  Eq.  (A12)  contribute  only  for  k  *^0,  since  they  are  independent  of  i;. 
However,  since  /4*"'s[y(t/2)]~'  is  known  by  definition,  we  are  only  interested  in  the  result  for 
A:  ^  I .  For  general  k,  we  can  make  use  of  the  relation 


/dV  I  _  y  •  ‘ 

\An)  i{fi  +  \)-p,  ^  (fin  +i)-F.)*’‘' 


and,  upon  taking  the  limit  »  0,  we  get  the  final  result 

1  I 


(A  1 3) 


(A  14) 


It  is  straightforward  to  prove  that  the  are  real  by  considering  the  symmetry  of  Eq.  (A  14)  with 
respect  to  the  substitution  n  -*  — r,  and  therefore  the  jj  are  real  also  [see  Eq.  (A17)]. 

Inserting  Eq.  (A14)  into  (A8),  rearranging  sums  and  defining  m  —k  —  \  gives 


r(p.)  (1/2  -  p,y  .?oC  ««  X  -  pJ ' 

This  relation  can  be  simplified  by  using  the  binomial  theorem  to  rewrite  the  factor 

ierXuT^)'=[-^r' 


The  result  is 
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.t*.  T'(P.)  l+4pi 


jSin(2A)C/  >  1). 


(A  15) 


(A16) 


(A  17) 


where  d,  is  defined  in  Eq.  (SO),  and  the  sum  has  been  restricted  to  positive  integers  only  (n  ^  I) 
by  use  of  the  symmetry  relation  p_,  =  -p,. 

APPENDIX  2 

Calculation  of  Expansion  Coefficients  for  Uniform  Solar  Disk 
For  the  case  of  a  uniform  solar  intensity  distribution,  Eq.  (49)  takes  the  form 


^  Z  77^  2  sin(2A) 

*  -  I  JoiPk)  I  +  '■P* 


(A  18) 


for  j  >  1,  with  70  =  [4//(l/2)]”'.  Writing  /ofP*)  =  (-  l)*l>fo(P*)l  and  defining  the  coefficient 

Pk  sin(2y0*) 


=  2 


l-^•(p*)l  1  +4pr 


(A  19) 


we  can  rewrite  E  ^.  (A  18)  as 


7>  =  2  I  (-l)V*. 
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(A20) 
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Generation  of  Trace  Constituent  Database 
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MODELING  OF  MILLIMETER-WAVE  EMISSIONS  OF  THE 
MIDDLE  ATMOSPHERE  OBSERVED  FROM  THE  GROUND  AND 
DEVELOPMENT  OF  CONSTITUENT  RETRIEVAL  ALGORITHMS 
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REMOTE  SENSING  AT  MICROWAVE  FREQUENCIES 


Satellite  and  ground-based  observations  at  microwave  frequencies  have  proven  to  be  the 
method  of  choice  for  determining  abundances  of  important  species  in  the  stratosphere  and 
mesosphere.  These  observations,  combined  with  a  solid  and  comprehensive  theoretical 
framework,  are  making  an  important  contribution  to  an  understanding  of  the  middle  atmosphere. 
CPI  has  been  providing  theoretical  support  to  NRL  in  this  area  through  the  development  of 
forward  model  calculations  and  retrieval  method  applications. 

In  the  summer  of  1989  Ms.  Donna  Kriebel  of  CPI’s  Annandale  office  began  working  with 
Dr.  Richard  Bevilacqua  of  NRL  on  the  early  phases  of  this  effort.  During  this  period  the  first 
versions  of  the  META  (Millimeter  Emission  Temperature  of  the  Atmosphere)  code  were 
produced.  At  the  same  time  a  congenial,  menu-driven  setup  code  was  developed  to  facilitate 
selection  and  modification  of  required  input  parameters. 

META  is  an  algorithm  which  generates  synthetic  emission  spectra,  using  radiative  transfer 
theory,  to  simulate  ground-based  microwave  measurements.  The  program  avoids  some 
approximations  of  earlier  forward  model  calculations,  and  has  an  added  advantage  in  that 
flexibility  has  been  a  major  consideration  in  its  design.  For  example,  it  can  be  used  in 
conjunction  with  appropriate  driver  routines  to  model  various  observing  modes,  such  as  beam¬ 
switching  or  single-beam  in  single  or  double  sideband  configurations.  Provision  is  also  made  for 
determining  background  effects.  Inputs  include  temperature,  pressure,  and  mixing  ratio  profiles, 
in  addition  to  the  characteristics  of  the  "measuring  system".  A  spectral  catalog  compiled  by  the 
Jet  Propulsion  Laboratory  is  consulted  for  selection  of  contributing  lines  according  to  user- 
specified  criteria.  The  simulated  spectra  so  generated  can  then  be  used  in  the  evaluation  and 
application  of  various  retrieval  schemes. 

Improvements  to  the  initial  version  of  META  include  a  driver  program  which  permits 
generation  of  a  number  of  spectra  using  different  mixing  ratio  profiles,  or  different  temperature- 
pressure  profiles,  calling  only  those  subroutines  which  are  necessary  to  take  into  account  the  new 
inputs.  Improvement  of  the  partition  function  entailed  expansion  to  include  vibrational  and 
electronic  parts,  whereas  previously  it  was  comprised  of  only  a  rotational  part.  In  some 
molecules,  such  as  ozone  and  CIO,  these  vibrational  and/or  electronic  states  assume  some 
importance.  In  addition,  a  new  pressure-broadening  half-width  was  incorporated  to  describe  the 
water  vapor  lineshape  more  accurately. 

In  the  ongoing  process  of  code  refinement,  comparisons  were  performed  in  which  the 
effects  of  varying  altitude  resolution,  varying  criteria  for  individual  line  selection,  and  inclusion 
of  the  troposphere  were  studied,  as  well  as  differences  between  species. 

Early  calculations  have  supplied  useful  comparisons  with  independent  results  of 
investigators  such  as  Dr.  Brian  Connor  of  NASA.  The  program  has  been  routinely  applied  to 
emissions  from  H20,  03,  and  CIO,  and  is  capable  of  modeling  spectra  from  a  wide  variety  of 
other  species 
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